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Abstract 

We perform a first principle semiclassical quantisation of the general finite-gap solution 
to the equations of a string moving on R x 5 3 . The derivation is only formal as we do not 
regularise divergent sums over stability angles. Moreover, with regards to the AdS/CFT 
correspondence the result is incomplete as the fluctuations orthogonal to this subspace in 
AdSz, x S 5 are not taken into account. Nevertheless, the calculation serves the purpose of 
understanding how the moduli of the algebraic curve gets quantised semiclassically, purely 
from the point of view of finite-gap integration and with no input from the gauge theory 
side. Our result is expressed in a very compact and simple formula which encodes the infinite 
sum over stability angles in a succinct way and reproduces exactly what one expects from 
knowledge of the dual gauge theory. Namely, at tree level the filling fractions of the algebraic 
curve get quantised in large integer multiples of H = 1/vA. At 1-loop order the filling 
fractions receive Maslov index corrections of \h and all the singular points of the spectral 
curve become filled with small half-integer multiples of h. For the subsector in question this 
is in agreement with the previously obtained results for the semiclassical energy spectrum 
of the string using the method proposed in |hep-th/0703191[ 

Along the way we derive the complete hierarchy of commuting flows for the string in the 
Mx5 3 subsector which are generated by the Taylor coefficients of the quasi-momentum p(x) 
through Hamilton's equation. Moreover, we also derive a very general and simple formula for 
the stability angles around a generic finite-gap solution which may be used in the study of 
stability properties of solutions in the MxS 3 subsector. We also stress the issue of quantum 
operator orderings and whether or not a given ordering preserves integrability since this 
problem already crops up at 1-loop in the form of the subprincipal symbol. 
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Introduction 



The method of semiclassical quantisation in field theory has been extensively developed by 
many authors in the 70' s using different approaches [H El El IU El El E] (see also the books 
P, E] for a more or less complete survey and list of references) . The aim of all these methods 
is to give a quantum mechanical meaning to extended classical solutions of the field equations 
which already classically exhibit particle like properties. The role played by such non-trivial 
classical solutions in the leading order quantisation of any field theory is evident from the 
path integral which is dominated by classical solutions in the h — > limit. It follows then that 
the applicability of semiclassical methods crucially relies on an explicit knowledge of classical 
solutions. Yet for a generic field theory, very little can be said about explicit solutions to the 
field equations and in most general solution does not exist. When the field theory 

is classically integrable however, essentially everything is known about the classical theory 
and the most general solution can be constructed explicitly in terms of standard functions 
and finitely many algebraic operations. In this case the complete semiclassical spectrum of 
the theory can then be obtained by applying the methods of semiclassical quantisation to 
the general solution. 

It is now very well established that the Metsaev-Tseytlin action [TU] describing super- 
strings on AdSs x S 5 is classically integrable [TT] . in the sense that the theory possesses an 
infinite number of integrals of motion. This fact has been thoroughly exploited in the litera- 
ture p2] to completely classify the full set of classical solutions on AdS$ x S* 5 by assigning to 
every solution a finite genus algebraic curve which encodes its integrals of motion Ji, . . . , I n . 
However, the algebraic curve is not enough to uniquely specify the solution, which can be 
seen as follows. Since a given solution carries only finitely many non-zero integrals of motion 
Ii, . . . , I n it will be invariant under the action of all the other integrals of the theory. More- 
over, the solution breaks all the symmetries generated by Ii, . . . , I n and the action of these 
integrals on the solution will generate new solutions with the same integrals. Indeed, in the 
theory of finite-gap integration [331 EH ESI EH1 EH EH EH] , (finite-gap) solutions are shown to 
be in one-to-one correspondence with sets of algebro-geometric data which essentially consist 
of a finite genus algebraic curve equipped with a finite set of points called a divisor. The 
action of the moduli Ii,...,I n of a solution on the solution itself will act non-trivially on 
the divisor, thereby generating a new solution with different divisor. The divisor therefore 
encodes the different zero-modes of a given solution. 

The treatment of the zero-modes is an important part of any approach to semiclassical 
quantisation [HE1IE1E]- Indeed, if a classical solution has zero-modes then a naive semiclas- 
sical quantisation of the solution will fail. Consider a solution <f) c \ of a field equation derived 
from an action S [(/>], i.e. S'[<p c \] = 0, where ' denotes 5 /5(f). If v denotes an infinitesimal 
symmetry of the equations of motion, i.e. v(S' [</)]) = S"[4>](v(j)), and suppose that (p c \ is not 
invariant under the symmetry then it follows immediately that (i>0 c i) 7^ is in the kernel of 
the operator S"'[0 c i] which is therefore not invertible and so the propagator of the theory in 
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the background (f) c \ cannot be defined. The standard way around this difficulty is to treat 
the zero-mode directions separately using the method of 'collective coordinates'. In short, 
collective coordinates parametrise the zero-mode directions, namely the flat directions in 
field space, along which the wave function will tend to spread out in the form of a plane 
wave as a result of which the quantum counterpart of the solution c i will acquire dynamics 
along the collective coordinates. Generally one has to perform a change of variables in field 
space to include the collective coordinates among the set of field variables and this can often 
only be done implicitly. A nice feature of the finite-gap construction is that it naturally lends 
itself to the separation of zero-modes since the divisor, which plays the role of the collective 
coordinates, already appears explicitly in the finite-gap solution - no change of variables was 
required. 

The divisor of a finite-gap solution therefore plays a central role in determining its semi- 
classical spectrum. But although the algebraic curve is known in full generality for the 
AdSz, x S 5 superstring, the divisor has only been identified so far in the subsector K x S 3 
for which the explicit reconstruction of finite-gap solutions from the algebro-geometric data 
has been studied [T3], [HI [15]. The method of semiclassical quantisation as stated above 
can therefore only be applied directly in the subsector M x S 3 . In this paper we perform 
such a semiclassical analysis of bosonic string theory on R x S 3 from first principles. We 
do not attempt to include the fluctuations in the directions transverse to the subspace 
K x S 3 C AdS§ x S 5 for clarity and because we believe that the method presented here 
should carry over with few alterations to the full case of superstrings on AdS§ x S 5 once the 
divisor is known. The calculation therefore serves as a toy model for understanding from 
the finite-gap perspective the origin of the discretisation of the algebraic curve when leading 
order semiclassical corrections are included. Nevertheless, our result agrees for fluctuations 
within the WL x S 3 subsector with the semiclassical result sj of [16l [TTl [18] . 

In the remainder of the introduction we start by recalling the method of semiclassical 
quantisation a la Dashen, Hasslacher and Neveu [H El [3] when applied to the specific example 
of the breather solution in Sine-Gordon theory. We reformulate everything in a language 
that we hope will facilitate the conceptual understanding of the method in the finite-gap 
setting and in the last part of the introduction we give a sketch of the ideas developed in 
the paper. 



0.1 Semiclassical Sine-Gordon breathers 

Consider the example of the boosted breather solution in Sine-Gordon theory [21 El E] 

»,.(«,«) = ^tan- { (( "" /2j > 2 - "'y" - TO > /(1 - " 2 >*,' 1 . (0.1) 
y/X \cosh[((rm/27r) 2 -l)5(27r/r)(x-^)/(l-t; 2 )l]J 

1 See also [HI HOI HI] for earlier work on obtaining the fluctuation energies from the spectral curve. In 
particular [21 where the one- loop energy shift was computed in the Landau-Lifshitz model. 
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This is really a two parameter family of solutions parametrised by their proper period r 
and their velocity v, or equivalently by their energy E and momentum p. To compute the 
(possibly continuous) spectrum of the corresponding quantum states it is always simpler at 
first to put the system in a very large but finite box of length L by identifying x ~ x + L so 
as to make the spectrum discrete, and then take the infinite volume limit L — ► oo at the end. 
In this closed-loop world the breather solution (10. ip is periodic in t of period T provided r 
and v satisfy T = /t/(1 — v 2 )^ = mL/v with l,m 6 N. 

If we were quantising the kink, we could move to its rest frame in which it is static and 
study small fluctuations in terms of eigenfrequencies. However, the breather is a little more 
complicated since it is time dependent in its rest frame, and because time dependent solutions 
are not point-like in field space, we need a way to characterise perturbations of the orbit as 
a whole. As we will describe in appendix [Bj this is done by considering the perturbation of 
a specific point on the orbit, evolving that perturbation under the equations of motion for 
roughly the period of the underlying solution, and comparing the final perturbation with 
the original one. If the perturbation is stable then it will have merely rotated and the angle 
of rotation is called the stability angle. If instead the perturbation is unstable it will have 
grown exponentially in magnitude, which corresponds to the case of a complex stability 
angle. Finally, if the perturbation comes back exactly to itself, this means it describes a 
nearby periodic solution, and in general zero stability angles correspond to symmetries. In 
the case of the Sine-Gordon breather we therefore need to look for generic nearby solutions 
4>(x,t) = <j) TiV (x,t) + 5(f>. This perturbed solution won't be periodic in general, yet because 
the linearised equation 



is invariant under time translation by T we can always write its solution as a superposition 
of eigenfunctions of time translation 8<f>(x,t + T) = e~ lu 5<f)(x,t), where v are their stability 
angles. Notice that the Sine-Gordon equation is invariant under arbitrary space and time 
translations, but the breather solution (f> TiV is not. As a result, d(f> TtV /dx and d<j) TtV /dt are 
both zero-modes, i.e. perturbations with zero stability angles. In general, any symmetry of 
the action that is not a symmetry of the classical solution will give rise to a zero-mode. 

The task of finding nearby solutions to the breather is greatly facilitated by the fact that 
the Sine-Gordon equation is integrable, since we can use the Backhand transform to get new 
solutions from known solutions. In particular we can perturb our breather by adding a little 
breather of small amplitude on top of it (Figured]). Studying double breather solutions in the 




(0.2) 



Figure 1: Perturbing the breather by another small breather using the Backlund transform 
limit where the small breather has vanishingly small amplitude corresponds to a linearised 

4 



study of the Sine-Gordon equation around the breather solution. So integrability gives us a 
convenient way of writing down explicit solutions to the linearised equation (10. 2p from which 
the stability angles of the breather may be read off. 

Identifying the space of classical solutions with phase-space, for each r, v (or equivalently 
E, p) the breather solution (10 .ip is just a specific point in phase space. However, the existence 
of two zero-modes d(p T}V /dx and d(p T}V /dt for the breather solution indicates that it really 
belongs to a two parameter family of solutions with the same integrals of motion E,p. These 
are the space and time translated breather solutions 



Since all the other stability angles of the breather are real, when we include first order quan- 
tum corrections the wavefunction will want to localise around not one breather, but around 
the whole two parameter family (I0.3P of breathers by spreading along the flat directions, 
namely the xq and to directions. Along these directions the wavefunction will therefore be 
a plane wave, but since the t -direction is closed by periodicity of the breather solution the 
plane wave along it must have an integer number of peaks and troughs. In other words the 
change of phase of the wavefunction around this closed direction will have to be an integer 
multiple n of 2tt. Along all the other non-zero stability angle directions the wavefunction 
will decay rapidly and, intuitively, for states with higher excitation number rii it will extend 
further in these directions. The correct quantisation conditions encoding the semiclassical 
energy spectrum of the wavefunction localised around the family of breather solutions was 
first derived by Dashen, Hasslacher and Neveu [1] and can be expressed as follows. If we 
define the 'action' of the breather solution as 



Although the derivation of this formula is very complicated, it intuitively makes a lot of sense. 
In general the phase of the wavefunction in the semiclassial approximation is an action of the 
form ( I0.4ap so the first term on the right hand side of ( 10.4b|) can be seen to comes from the 
single-valuedness of the wavefunction along the compact ^-direction whereas the correction 
from the sum over stability angles is related to the small fluctuations transverse to the t 
and Xq directions. 

For the purpose of drawing the analogy between the Sine-Gordon breather case here and 
that of finite-gap strings discussed later it will be convenient to think of the conditions ( 10.41) 
in more geometric terms in phase-space as follows. Since the breather in (10.31) with xq = 



4> TtV (x + x ,t + t ). 



(0.3) 




(0.4a) 



then the DHN quantisation conditions read 




(0.4b) 
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is periodic, it can be thought of as a closed orbit on the level set of fixed E,p. The 
direction along the orbit, parametrised by t , corresponds to the zero-mode d<j) TtV /dt of the 
breather. But since it has another zero-mode, namely d(p TiV /dx, this orbit really belongs to 
a continuous family of periodic orbits, parametrised by xq, all contained in He, p - However, 
because we are working in a periodically identified finite box, this two parameter (xo,to) 
family of breathers is in fact a torus T Ep lying within And since all the other stability 

angles of the breather are non-zero, this means that T 2 E p is isolated on the level set Ee iP 
in the sense that it does not belong to a larger continuous family of periodic orbits within 
T,e, p - Yet if we leave the level set £e, p , one can show that in a neighbourhood of S^ iP the 
torus T E persists, namely it belongs to a two parameter family of torii parametrised by 
E,p. This is the content of the 'cylinder theorem', illustrated in Figure [2] for the case of 
a solution with a single zero-mode, so that its zero-mode family in the level set H~ l (E) is 
just a circle S E that belongs to a cylinder S E x [E — e, E + e] . Looking back at the most 




Figure 2: Cylinder theorem: a periodic solution 7^ on the energy level H 1 (E) is contained 
in a one parameter family of periodic solutions of varying energy in the range [E — e, E + e] . 

general breather solution (10.31) it contains four independent parameters: the two parameters 
£0, to are parameters along the torus T E whereas E, p parameterise the family of torii of the 
cylinder theorem. Now the effect of the quantisation condition (10.41) is to pick out a discrete 
set of breathers from this 'cylinder' of breathers ( 10 .11) . the energy and momentum of which 
approximate to order 0(h) the semiclassical energy spectrum of the quantum states localised 
around the breather solution. For instance, when applied to the Sine-Gordon breather the 
quantisation conditions (10.41) yield the following semiclassical spectrum [2] 

Ek,n = (pl + M^, Pk = ^, 

where M n = sin ^ and 7' = ^ (l — 87T X m2 ) \ and in the infinite volume limit L —> 00 
the momentum becomes continuous as expected. 

0.2 Sketch of semiclassical finite-gap strings 

We would like to apply a similar kind of reasoning to the case of superstring theory on 
AdS^ x S* 5 . However, since this formalism requires the knowledge of explicit solutions we 
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will restrict attention to bosonic string theory on R x S 3 for which the general finite-gap 
solution to the equations of motion is known [121 fH| - In conformal static gauge the string 
is given by an embedding g(a,r) G SU(2) of the worldsheet into SU(2), and if we define 
the corresponding Lie algebra current j = —g~ 1 dg G su(2) then the equations of motion and 
Virasoro constraints take the following form 

d*j = 0, dj- 3 Aj = 0, ^trf ± = -K 2 . (0.5) 

As is well know, the equations of motion are integrable and can be rewritten in the form of 
a zero-curvature equation dJ(x) — J{x) A J{x) = 0. In this form one can make use of the 
powerful methods of finite-gap integration to construct, at least abstractly the general finite- 
gap solution to the equations of motion. In fact, it is possible to incorporate the Virasoro and 
static gauge constraints into the constructions [131 E] so as to get only physical motions 
of the string. The general finite-gap solution is constructed from the following piece of 
algebro- geometric data: 

• An algebraic curve p2] of genus g. 

• A set of g + 1 points [13] on this curve. 

Essentially, by the Riemann-Roch theorem there is an injective map from this algebro- 
geometric data into the space of solutions to (10.51) . The idea of finite-gap integration is 
illustrated in Figure [3j Every finite-gap solution to (10.51) is in one-to-one correspondence 




Figure 3: Idea of finite-gap integration. 

with an algebraic curve (of genus three in Figure [H]) equipped with a set of marked points 
(four of them in Figure [3]). The algebraic curve encodes the integrals of motion of the 
solution, and these points encode the dynamics. Their exact motion on the algebraic curve 
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is very complex, but what we find is that if we map the algebraic curve to its (generalised) 
Jacobian, a (g + l)-torus, via the (generalised) Abel map then the motion in a and r becomes 
extremely simple, namely it linearises. The motion of the string on this (g + l)-torus is like 
that of an infinitely rigid string wrapping one cycle of the torus and moving linearly in time 
along another direction. 

An alternative way of picturing what a finite-gap solution looks like that will be useful 
later is as follows. As we just saw, the dynamics of a finite-gap solution corresponds to 
linear motion on a (g + l)-torus, which is very reminiscent of a finite-dimensional integrable 
system. In fact one can view the Jacobian as the Liouville torus of a (2g + 2)-dimensional 
dynamical system. The base space C of this (2g + 2)-dimensional system is the moduli space 

(0) 



Figure 4: The algebro-geometric data as a (2g + 2)-dimensional phase-space. 

of the algebraic curve parametrised by the filling fractions {Si = zdp} 9 ^. But if the 
algebro-geometric data is to be thought of as a finite-dimensional phase-space it must be 
equipped with a natural symplectic structure. This can be obtained as follows: the finite-gap 
solution maps this algebro-geometric data to the space of solutions to ( 10. 5ft . see Figure HI 
Identifying the space of solutions to the equations of motion with the phase-space V°°, the 
solutions to ( 10. 5ft which also satisfy Virasoro and static gauge define a second class constraint 
surface V v C This is equipped with a Dirac bracket induced by the Poisson bracket 
on the su(2) current appropriately regularised a la Maillet [221 E]. If one then pulls back 
this Dirac bracket to the algebro-geometric data using the finite-gap solution we obtain a 
'natural' symplectic structure on the algebro-geometric data which can be concisely written 
as (see [H] for details) 

3+1 

u> = dSj A dfj. 
i=i 

The upshot of this is that the filling fractions are precisely the action variables of the finite- 
gap string. They are the analogues of the period r and velocity v (or energy E and momen- 
tum p) of the generic breather ( 10.31) which defined a four parameter family of solutions. A 
finite-gap solution defines a whole (2g + 2) parameter family of solutions parametrised by 
the algebro-geometric data and can be written as follows 



g = g(J2 N t " u ^ + D \ s ) 



where are a set of g + 1 independent times (defined in section [2]), S is the vector of 
action variables which plays the role of the parameters (r, v) here and D £ <C 9+l is the exact 
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analogue of the initial coordinates of the breather (xQ,t ). We therefore expect a finite-gap 
solution constructed from a curve of genus g to have g+1 zero-modes corresponding to the 
g + 1 components of the vector D. 

In view of applying a semiclassical quantisation formula like the one in (10.41) we must first 
determine all the stability angles of a given finite-gap solution. So just as in the case of the 
Sine-Gordon breather, we would like to study perturbations of finite-gap solutions described 
above. Once again integrability will play a prominent role in solving the linearised equations. 
In fact, finding solutions to the linearised problem is very simple now that we have already 
fully exploited integrability to construct the most general finite-gap solution. A perturbation 
of a given finite-gap solution will simply be another 'nearby' finite-gap solution. Recall 
[j~3l [T4"] that in the SU (2) sector the algebraic curve is hyperelliptic and can be represented 
by a set of g + 1 cuts in the complex plane. How can one describe perturbations of the 
g-gap solution corresponding to this curve? Playing the same game as for the Sine-Gordon 
breather where we used integrability to add another little breather on it, here we can just 
take a solution corresponding to a curve of genus one higher, but make the extra filling 
fraction very small, which corresponds to making the cut very small, see Figure El There 



is an obvious analogy here between breathers in Sine-Gordon and cuts in bosonic strings 
on M x S 3 as one can think of a finite-gap solution as a multi-breather solution consisting 
of finitely many breathers. Cuts with small filling fractions are analogous to breathers of 
small amplitude as both describe perturbations. If we define the a^-cycle {i = 1, . . . ,g) as 
usual to encircle the i th cut counterclockwise (on the upper sheet) then a perturbation of 
this kind clearly corresponds to pinching an a-cycle of the algebraic curve. So we want to 
take the difference between the solution before pinching an a-cycle and the solution after 
pinching the a-cycle; this will give us a perturbation of the latter and we can then analyse 
its periodicity properties to extract the corresponding stability angles. Notice however that 
any given perturbation of a finite-gap string will have one stability angle defined for each 
cycle on the Jacobian, or equivalently for each macroscopic cut. 

The semi-classical spectrum can be obtained by performing a WKB analysis of the wave- 
function that will localise around the zero-mode directions of the solution, which in the case 
of the finite-gap string is the Jacobian. Again, the leading term will describe how many full 
waves fit on the compact Jacobian, and the infinite sum corresponds to small fluctuations 




: 



Figure 5: Perturbation of a finite-gap solution. 
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transversal to the Jacobian. The result of such an analysis that will be sketched in section 
11.21 are the following set of Bohr-Sommerfeld equations, the correct form of which involves 



Here \ii — 2 is the Maslov index of the ^4/-cycle (I = 1, . . . , g+1) in the generalised Jacobian 
J(S, oo^). Note that (10.61) is only valid in the harmonic oscillator approximation Nj ^> n a 
where the perturbations are much smaller than the background filling fractions. So the 
expression ( 10. 6p really contains two different orders, namely the tree level and 1-loop level 
of order 0(1) and 0(h) respectively (after multiplying ( 10.61) throughout by h). At tree level 
(10.61) simply expresses the fact that the filling fractions are quantised in integer multiples 
of h, i.e. Sj = Njh, which is a straightforward consequence of the fact that the Si are the 
action variables as was shown in p3]. The non-trivial content of ( 10. 6ft is the 1-loop correction 
which includes firstly the Maslov index correction ^-h and secondly the infinite sum over 
stability angles. 

Obtaining the energy spectrum from (I0.6P is relatively straightforward since for a system 
to be semiclassically integrable requires that [Si, Sj] = 0(h 3 ) and so the energy eigenvalues 
are given to leading order in h simply by evaluating the classical energy E C \[S\, . . . , iS^-j-i] on 
the eigenvalues of the action variables ( 10.61) . As we show in section H~T1 this can be expanded 
to order O(h), expressing the result as a sum of the tree level term E c \[Nih, . . . , N g+ ih] and 
the 1-loop correction involving the sum over stability angles [T6l [T71 [T8] . But moreover, in 
section 14.11 we also show, using the result of section 12.31 that the sum of the tree level term 
and 1-loop correction term can be succinctly rewritten in a compact form that captures the 
complete result at 1-loop in a unified way. Indeed, we show that the energy spectrum can 
be formally obtained by evaluating the classical energy of an infinite-gap solution with all 
its infinite filling fractions quantised to half-integer multiples of h, namely 



This result is to be interpreted as a limit of expressions where a finite but arbitrary number 
of first entries are of order 0(1) corresponding to the tree level order and the remaining 
infinite number of entries encode the stability angle contribution to the 1-loop corrections of 
order O(h), see (I4.5p . 

The paper is organised as follows: in section [1] we review some basic features of semiclas- 
sical quantisation for finite-dimensional systems. In particular we remind the reader how 
operator ordering enters in the semiclassical regime: in the language of pseudo-differential 
operators (appendix [A} the different operator orderings are encoded in the subprincipal 
symbol [31]. We also sketch the derivation of the Bohr-Sommerfeld quantisation conditions 

2 Here and in the remainder of the paper, in the string theory context we will always let H = -4=. 





(0.6) 
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[2H [221 EE] • In section [2] we look back at the general construction of finite-gap strings P2HH] 
and derive the whole hierarchy of commuting flows. That is, we show how the integrable 
equations of motion for the embedding of the string in K x S 3 are part of an infinite hierarchy 
of higher integrable equations corresponding to the infinite set of conserved charges of the 
string, as is usual in any integrable system. In section [3] we discuss perturbations of a generic 
finite-gap string through the pinching of a-cycles. This leads to a general formula for the 
non-zero stability angles of a generic finite-gap string. Using this result, in section H] we come 
back to the issue of semiclassical quantisation of finite-gap strings and apply the formalism 
of section [1] to obtain the semiclassical spectrum of the string. Some appendices elaborate 
on the discussion in each section. 



1 Semiclassical approximation generalities 

Consider a classical Hamiltonian system described by a 2n dimensional phase-space T*X 
with Hamiltonian H : T*X — > R. Given Belwe can consider the codimension one energy 
level set = H~ l (E) C T*X. Assume also that we have a desired quantisation of the 
system, that is, we have a self-adjoint operator H acting on L 2 (X) whose principal symbol is 
the classical Hamiltonian H. If H~ 1 ([E — e, E + e]) is compact then the eigenvalues of H in 
the range [E — e, E + e] will be discrete since the corresponding eigenf unctions are localised 
around this compact set. The goal of semi-classical quantisation is to obtain the spectrum 
of H in [E — e, E + e] to leading order in h. One approach is to describe the spectrum using 
what are known as trace formulae, the basic idea being to encode the spectrum in terms of 
a single function n(E) = ^°1 S(E — Ef) = tr 5(E — H) where Ej denote the eigenvalues 
of H and which can be rewritten as 

n(E) = Re^- dt tr e ^ E ~ H)t = Re^- / dt [^e"^. (1.1) 

7rhJ TrhJ J p-o. 7 

period t 

In the semiclassical limit S^Owe perform a stationary phase approximation of the integral 
on the right hand side in order to obtain a semiclassical estimate of the spectrum {Ej 1 } 
of H. The presence of the trace means that dominating contributions come from periodic 
orbits of the classical system. This is a general feature of semiclassical trace formulae which 
relate analytic data of the operator H (namely its eigenvalues) to geometric data of the 
corresponding classical Hamiltonian H (namely its periodic orbits). This is one advantage 
of trace formulae over other semiclassical quantisation methods in that they elucidate the 
relation between the semiclassical spectrum and the classical periodic orbits. 

On the downside however, despite the geometrical appeal of the path integral approach 
to semiclassical quantisation, it is hard to discuss the issues of operator ordering within this 
framework. Indeed, thinking in terms of phase-space path integrals, since everything in the 
integrand itself is classical, any information about quantum ordering is neatly tucked away 
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in the definition of the regular isat ion used in the phase-space path integral measure [dj]. 
The standard choice of discretisation of the path integral measure involves the mid-point 
prescription which corresponds to the Weyl-ordering prescription in the operator formalism. 
In particular the quantum Hamiltonian is the Weyl-ordered classical Hamiltonian, i.e. H = 
Op^(H). In order to deal with operator ordering issues, it is therefore more convenient to 
work directly with operators. 

A convenient operator formalism for discussing semi-classical quantisation involves pseudo- 
differential operators (referred to as \l/DOs for short). We refer to appendix |X] for a very 
brief introduction to \l/DOs and their relevance for treating semiclassical quantisation. The 
basic idea of this approach is to associate with any operator / not a single function on T*X, 
which cannot by itself encode all the information about operator ordering in /, but a family 
of functions fn, € C(T*X) called symbols. The leading function f is exactly the classical 
function corresponding to /, whereas all the sub leading functions encode the operator order- 
ing in /. So instead of working with operators one can work directly with their respective 
symbols. Moreover, in the semiclassical approximation one only needs to deal with the first 
two symbols of an operator, known as the principal symbol (i.e. the classical function) and 
the subprincipal symbol. We will turn to the formalism of ^DOs and the issue of operator 
ordering in an integrable system in section 11.11 In section 11.21 we will show how the Bohr- 
Sommerfeld quantisation conditions are modified by the presence of a subprincipal symbol 
which reflects a choice of ordering. 

But first, to get an intuitive idea of how operator ordering ambiguities arise even at 
the semiclassical level to affect the quantisation conditions, it is instructive to consider the 
simple example of the harmonic oscillator for which the leading order quantisation is exact. 
The classical harmonic oscillator Hamiltonian is H = £ — h hmuj 2 x 2 , and the action variable 

2m 2 ' 

of the closed path of energy E is given by 

I = — (f pdx = — . 
2vr J h=e uj 

By promoting the variables x,p to operators x,p there is only one reasonable choice of 
ordering in the Hamiltonian, namely the Weyl-ordered Hamiltonian H = ^— + |mu; 2 x 2 . 
The spectrum of such an operator is well known to be E n = (n + |) frxj,n G N so that 
the spectrum of the Weyl-ordered action variable I = -H is simply given by the standard 
Bohr-Sommerfeld quantisation condition, 

Spec (I) C {% + ij h, 

where the index of \ by which the spectrum is shifted from KL is known as the Maslov 
index in the context of Bohr-Sommerfeld quantisation. Now since we are given at the outset 
only the classical Hamiltonian, we could always choose to quantise it with a more perverse 
choice of ordering. For instance, if we rewrite the classical Hamiltonian as H = uaa* where 
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a = yY^jr (% + <fj^) and after promoting everything to operators request that in the quantum 
Hamiltonian the a sits to the right of the a) then we obtain the normal-ordered Hamiltonian 
:H: = ouHcra, where [a, a'] = 1. The corresponding normal-ordered action operator is given 
by : I : = ha) a whose spectrum is easily seen to consists of integer multiples of h, 

Spec (:/:) C ZH. 

We observe that the Maslov index is precisely cancelled by the shift from Weyl-ordering to 
normal-ordering. Even though in the case of the harmonic oscillator we know that the correct 
physical quantisation of H is the Weyl-ordered one H we would like to stress that in general 
the choice of operator ordering in the quantisation of the action or Hamiltonian may not 
be as obvious and their spectrum may observe a shift from the standard Bohr-Sommerfeld 
spectrum (Z + jj h, where ji G Z 4 is the Maslov index. 

1.1 Operator ordering issues 

As explained in appendix |XJ one can keep track of operator orderings in the language of 
pseudo-differential operators by retaining subleading terms beyond the principal symbol 
in the full Weyl symbol of an operator. In most applications of the theory of \l/DOs the 
quantities of interest are specified as \l/DOs at the outset so that their full Weyl symbol is 
known. In the present case however we start from a classical system specified by its phase- 
space (T*X, uo) and the set of classical observables of interest are iq, . . . , F n , H. Quantising 
this classical system requires an operator ordering prescription for obtaining operators from 
the corresponding classical observables. At the semiclassical level this boils down to the 
specification of an extra function, the subprincipal symbol, for each classical observable. 
Specifically, given a classical observable /o £ C(T*X), we construct 

/=o P r(/ +M), 

where the presence of the subprincipal symbol fi G C(T*X) reflects the operator ordering 
ambiguities already manifesting themselves at the semiclassical level. Every possible choice 
of a function fi G C(T*X) corresponds to a different prescription for the operator ordering in 
/ at order 0(h). The principal symbol fo = GqQ) is the corresponding classical observable. 

Recall the definition of an integrable system, which roughly speaking is one which pos- 
sesses the maximum possible number of independent integrals of motion. Specifically, a 
Hamiltonian system (T*X, H) is said to be classically integrable if there exists n functions 
F u . . . , F n G C(T*X) such that 

(1') dF\ A ... A dF n ^ almost everywhere, 
(2') {F i ,F j } = 0, Vi,j = l,...,n, 
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(3') H = H(F 1 ,...,F n ). 



Conditions (2') and (3') together imply that the Fj are in fact integrals of motion, XjjFi = 0. 
In other words, T*X admits a torus action with moment map 

F = {Ft,...,F n ) :T*X^R n . 

At regular values / of F, the level sets F~ l (f) define n-torii (in the compact case) and 

F 

foliate the phase-space T*X; namely T n T*X — > M. n . This foliation allows one to define 
canonical action-angle coordinates with the action variables {/j}™ =1 parametrising the base 
M. n and the conjugate angle variables {#i}™ =1 , each taking values in [0, 2ir}, parametrising the 
independent cycles of the torus T n . 

We will say that a \l/DO H is semiclassically integrable if there exists n \l/DOs F\, . . . , F n 
with principal symbols Fj = Oq (Fi) such that 

(1) dFi A ... A dF n ^ almost everywhere, 

(2) [Fi, Fj] = 0(h 3 ), Vz,j = l,...,n, 

(3) H = H(Fi, ...,F n ) + 0(H 2 ) for some function H. 

Notice that we only require commutativity modulo 0(h 3 ) in property (2); it guarantees 
in particular that the operator H(F±, . . . , F n ) in (3) is free of operator ordering ambiguities 
certainly up to 0(h 3 ), so that property (3) makes sense. Property (2) is to be contrasted with 
the definition of full quantum integrability which requires exact commutativity [Fi, Fj} = 0. 
Now since a^([Fi, Fj}) = -ih{F h Fj} (see appendix [A]) and a^(H) = H(F U . . . , F n ), it 
follows that the principal symbols Fi = ' (Fi) satisfy all three properties (l')-(3') above 
for a classically integrable system with Hamiltonian H = Gq(H). This means that any 
semiclassically integrable system exhibits at leading order the full geometric structure of the 
underlying classically integrable system given by its principal symbols. In particular, the 
level set Aj = F~ l (f) of the moment map F = (F±, . . . , F n ) : T*X —> M n is a Lagrangian 
n-torus and foliates phase-space T*X as we let / vary. 

But the notion of semiclassical integrability contains more information than that of its 
underlying classical integrable structure [3H[32]. Property (1) only contributes at leading 
order since it is a statement about the principal symbols Fj alone, whereas property (2) at 
0(h 2 ) yields an equation for the subprincipal symbols F? = cr^ b (Fj) of the Fj (see appendix 

ED 

= jf^FiM = {Fi,F°} + {F°,Fj}. (1.2) 

It is possible to interpret these equations geometrically so as to supplement the geometrical 
structure already laid out by the principal symbols with further geometrical data. For this 
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we define the subprincipal form k on Aj- by denning its action on the basis vectors X F . at 
any point of Af through [3_T] 

K(X Fi ) = -F*, i = l,...,n. (1.3) 

It then follows immediately from (11.21) that k is closed since 

dK(X Fi , X F .) = X Fi K{X F .) - X F .K(X Fi ) - k([X Fv X f .}) 

= -X Fi F° + X Fj F° - k(X {Fi , Fj} ) = -{F u F°} + {Fj, F*} = 0. 

Hence the operator ordering in the Fi can be accounted for at the semiclassical level by 
specifying a closed 1-form k on the Liouville n-torus Af. And in fact it is clear from (11.31) 
that every choice of a closed 1-form k G Q 1 (Af) corresponds to a different choice of operator 
ordering in the definition of the F^. 



1.2 Bohr-Sommerfeld conditions 

We are interested in the joint spectrum of the Fi up to 0(h) which requires solving the 
eigenvalue problem to that order 

(F l -/ l )^ = 0(n 2 ). (1.4) 

The Bohr-Sommerfeld conditions are conditions for the existence of a solution to these 
coupled pseudo-differential equations. Their rigourous derivation is rather involved but here 
we would just like to outline how the subprincipal symbol comes about in these conditions. 
To solve ( 11 Ah locally one considers a local patch V C Af on which n : T*X — > X is a 
diffeomorphism and uses the WKB ansat^l 

^WKB = e^-^°p + 0(h) 

on U = 7r(V) C X where the nature of p will be specified shortly. If we let : U T*X 
denote the 1-form c£0_i viewed as a map then equation (ll.4j) implies to leading order in h 
that [23] 

im = V C Af, (1.5) 

or i>d4-i — My 1 - a P ro P er ty of the tautological 1-form a, namely = we then 

have [23] 

dn\* v <p~i = a, (1.6) 

3 The WKB ansatz isn't actually restrictive since one can show that the space of solutions to (|1.4|) is one 
dimensional so that any other solution is proportional to the WKB solution. 
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in other words, 7r|y0_i is a local solution to the classical integrability condition uj = da = 
on Aj. If p is a half-densitjH on U C X then the subleading order of (jl.4p can be written 
invariantly as (Theorem 11.11 pl26 of [21]) 

(-iC XFi +F°) Ht^V) =0. 

Now provided the subprincipal symbols are real this equation implies on the one hand that 
n\yp is an invariant half-density on Aj, i.e. Cx F ^\vP — 0, and on the other hand that 

dTt\v<po = k, (1.7) 

which says that 7r|y0o is a local solution to the subleading integrability condition that dn = 
on A/. What one would like to do is patch up the local WKB solutions ipwKB defined on 
local neighbourhoods of Ay for which the projection tt : T*X — > X is a diffeomorphism. 

A f C T*X 

singularity 



Figure 6: Caustics of the Lagrangian submanifold Af 

However, one runs into problems at caustic points where tt is singular (see Figure E]). A 
way around this problem was proposed by Maslov which allows one to define a solution to 
( 11.41) which is localised and defined patchwise on Ay (near caustics one uses the "momentum" 
projection ir p of T*X onto a typical fibre of T*X instead of tt). The single valuedness of this 
global solution requires its phase to be an integer multiple of 2tt. The phase is essentially 
that of the local WKB solutions ipwKB introduced above but with additional Maslov index 
corrections (coming from the caustics) so that single valuedness conditions, known as the 
Bohr-Sommerfeld-Maslov conditions, read 

^hl a + ^L K = N ' + t T + 0ihl i = i -- n (L8) 

where 7^ is a basis of ifi(Ay,R) with Maslov indices /i 7i G Z 4 and integers Ni G Z. Note 
in particular the presence of the subprincipal form k which as we have argued is related to 

4 Since the product of two half-densities is a density of weight one there is a natural inner- product on half 
densities (pi, P2) = J M P1P2 which makes the completion into a Hilbert space. 
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operator ordering ambiguities in going from a classically integrable system to its quantum 
(or just semi classically) integrable counterpart. It has the effect of shifting the spectrum 
of the action variables similar to what happens in the case of the harmonic oscillator when 
we change quantisation, from Weyl to normal ordering say [27] . In the cases where all 
the operators are chosen to be Weyl ordered, in particular the Fi, we have k = and (11.91) 
reduces to the EBK quantisation conditions. In the remainder of the paper we shall make the 
assumption that the cohomology class [k] e if 1 (Ay) of the subprincipal form k vanishes. The 
reason for this assumption is that the result is simpler to express in this case and moreover it 
agrees with the results of [TBI C3 EE]. With this assumption, the Bohr-Sommerfeld-Maslov 
conditions simplify 



We stress that this assumption does not imply the choice of Weyl ordering since it only 
corresponds to setting the subprincipal symbol to zero, whereas Weyl ordering corresponds 
to setting all the lower order Weyl symbols to zero as well. 

Now the derivation of the Bohr-Sommerfeld-Maslov conditions (11.81) or (11.91) essentially 
consisted in quantising a Lagrangian n-torus Af by constructing a wave- function localised 
around it. However, even though the level set Af = i r ' _1 (/) is indeed a Lagrangian n-torus 
for almost every value of the integrals of motion fi, . . . , f n in an integrable system, there 
exists interesting level sets F~ 1 (f) in phase-space where this is not the case. This happens at 
the (measure zero) set of critical values of the map F = (Fi, . . . , F n ). Consider for instance 
the two-dimensional harmonic oscillator with different frequencies and Hamiltonian 



whose integrals of motion are given by Hi,H 2 . For non-zero values E\ , E 2 ^ of Hi , H 2 
the level sets H~ l (Ei,E 2 ) consists of two ellipses, in other words a Lagrangian 2-torus. 
However, if say E 2 = the level set H~ 1 (Ei,0) consists of just a single ellipse (Figure [TJ. 
The same thing is true when E\ — and at the point where E\ = E 2 = the level set 




(1.9) 



tt Pi , 1 2 2 i P2 1 1 2 2 rr,rr 

H = 17 + 0^1 + -W + 0^2 =Hi + H 2, 



(1.10) 



pi 




Figure 7: Periodic orbit with H 2 = of energy H = H\ = E. 
consists of just a single point. One can draw a picture of the phase-space in the region 
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where £ = {(E 1 ,E 2 ) : > 0, i — 1,2} which is foliated by 2-torii in the interior of S 
but with the fibres over the boundary dS \ {(0,0)} being ellipses and the fibre over the 
point (0, 0) being just a single point, see Figure El Note that the set of critical values 




Figure 8: The phase-space of the two-dimensional harmonic oscillator. 



88 is of measure zero. However, if we are interested in the semiclassical spectrum of the 
two-dimensional harmonic oscillator in the region near 88 then a modification of the Bohr- 
Sommerfeld-Maslov quantisation conditions (11.91) is required so that it applies to isotropic 
p-torii which are the level sets of a limited number p < n of integrals of motion iq, . . . , F p . 

It was pointed out by Voros ES] that the Bohr-Sommerfeld conditions (11.91) for the 
apparently more restrictive case of an integrable system may be used to obtain the Bohr- 
Sommerfeld conditions in all other intermediate cases, namely the partially integrable one 
(with p < n integrals of motion) and even the non-degenerate case p = 1 (where H is the 
only integral). If the system has p independent observables F = (iq, . . . , F p ) in involution 
(with H = H(F)), then on each codimension p level set = F^ x (f) the system has a p- 
torus A/ C £/ generated by the vector fields Xp v Each of these p-torii is surrounded by an 
n-torus of the linearised system to which the Bohr-Sommerfeld-Maslov conditions ( 11. 9ft may 
be applied. This results in a set of Bohr-Sommerfeld conditions for the cycles on the p-torus 
which include stability angles for the small fluctuations in the directions transverse to this p- 
torus. The derivation of these Bohr-Sommerfeld conditions from those in the integrable case 
( 11.91) are a bit lengthy but the derivation in the more general case 1 < p < n is conceptually 
the same as the p — 1 case [25j [26]. For completeness and to explain the appearance of 
the stability angles (which are related to the eigenvalues of the Poincare map) we repeat 
the details of the derivation of [23 EE] in appendix [B] The result (IB. 1 [) is the following 
quantisation condition for the isolated orbit 7 



h + o(h 2 



1.11: 
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where N G Z and n a G N. Since the periodic orbit 7 C Eg in fact belongs to a continuous 1- 
parameter family j(E) of periodic orbits parametrised by the energy E ('cylinder theorem'), 
what the condition (11. lip does is pick out a discrete set of periodic orbits j(E^), in a 
neighbourhood of Eg, whose energies i^j 1 approximate eigenvalues of H to leading order 
in h, see Figure [9j The condition depends on the stability angles v a G R (defined via 



Figure 9: Bohr-Sommerfeld semi-classical spectrum: the discrete set of periodic orbits r ){E^ j 
shown in blue have energies E^ a approximating the eigenvalues of H to 0(h 2 ). 



the Poincare map, see appendix [B]) of the stable isolated orbit 7 and is valid only in the 
approximation where < <§; |A^| which is required for the linear approximation (used in 
deriving these condition) to hold. 

The more general case of a system which has p independent observables Fi,...,F p in 
involution (with H = H(F\, . . . , F p )), where p lies in the range 1 < p < n is a straightforward 
generalisation. In this case we get a set of p quantisation conditions, one for each cycle 
jk, k — 1, . . . ,p on the p-torus 



/ a 



27T (N k + ^ 



4 J 



E 

a=p+l 



1 » 

2 A 



h + o(h 2 ). 



(1.12) 



This time there are p conditions on the p parameters fi, ■ ■ ■ , f p of the codimension p level 
sets E f = F- x (/). 

To illustrate the use of the modified Bohr-Sommerfeld conditions (11. lip for an isolated 
orbit let use go back to the case of the two-dimensional harmonic oscillators (11.101) . This 
system is obviously integrable and the exact spectrum of H is given by 



ni + - I hwi 



f - ) ^2- 



However, suppose for the sake of argument that we can only solve classically for the Hamilto- 
nian H\ and wish to obtain the spectrum of H — H\ + Hi by perturbation as describe above. 
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Then consider a particular motion of the Hamiltonian Hi of total energy Hi = E, through 
the point (pi,Xi,p 2 , x 2 ) = (p , 0, 0, 0) say, see Figured This defines a 1-parameter family of 
periodic orbits parametrised by their energy H = Hi = E. It is clear that the (p2, x 2 )-p\a,ne 
gives a Poincare section of the orbit through the point (po, 0, 0, 0) since all orbits of Hi have 
the same period T\ = The prescription for determining the stability angles of this orbit 
is to consider small perturbations around it within the same energy level H = E. If the 
periods of the two harmonic oscillators are different, Ti ^ T 2 , then after a length of time Ti, 
the motion in the (p 2 , a^-plane does not close and there is a deficit angle of v — lo% ■ T\ , see 
Figure [101 The tower of energy levels corresponding to the periodic motion in Figure H is 

PI P2 




Figure 10: Perturbed trajectory of energy H — Hi + H2 — E. 
therefore given by the Bohr-Sommerfeld condition (11. lip which in this case reads 



h 



rii + 



n 2 + 



2 2vr 



h + o{h 2 ) 



and hence E ni>n2 — U\-I\ = (rii + |) t\wi+ (n 2 + |) huj 2 +0(h 2 ) so that the Bohr-Sommerfeld 
condition is actually exact to first order in h on the harmonic oscillator. 



2 The string hierarchy 

In the general theory of finite-gap integration j33j EH [351 ES E3 EEJ [39] the reconstruction 
of a solution requires an algebraic curve E, which specifies the integrals of motion, as well 
as a divisor (i.e. a finite set of points) 70 on S which specifies the initial conditions for the 
dynamics. However, since the system is integrable it possesses an infinite number of integrals 
of motion H{, each one generating a different Hamiltonian flow on phase-space in the usual 
sense through the Hamilton equations dtj = {Hi,f}. The dynamics of the divisor 7(tj) 
(with initial condition 7(0) = 70) on E with respect to some time ti is then determined by 
the corresponding Hamiltonian Hi. In fact, in this setup there is a natural correspondence 
between Hamiltonian flows and meromorphic differentials on the algebraic curve E. This 
is well known for instance in the case of the worldsheet coordinate a which couples to the 
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quasi-momentum dp. Indeed, it was shown in [T3] that the worldsheet coordinates (er, r) 
enter the finite-gap solution only through the meromorphic differential 

dQ — (adp + rdq) , 

so that the coordinates a and r are said to "couple" respectively to the quasi-momentum 
dp and the quasi-energy dq. The aim of this section is to similarly identify the dynamics 
corresponding to all the higher conserved charges within the finite-gap language in terms of 
meromorphic differentials on the underlying curve E. 

The reason for doing this is the following. As we have already mentioned in the introduc- 
tion and will recall again in section HI a finite-gap solution can be understood as an injective 
map from a finite dimensional phase-space to the full infinite dimensional phase-space of the 
theory [HJ [35] . The Liouville torus of the finite dimensional phase-space in question is the 
(generalised) Jacobian J(£, cxd ± ) of the algebraic curve £ which is a (g + l)-torus. We will 
show that the divisor moves linearly on J(£, oo ) with respect to all the higher flows. But 
since J(E, oo*) is g+ 1 dimensional one can use g + 1 independent such flows to parameterise 
it. This will give a nice coordinate system on the Jacobian which will be useful when we come 
to consider perturbations of this (g + l)-torus in section [3] for computing stability angles. 
In particular, the angle variables (fj which will couple to the quasi-actions dq^ defined in 
section [2~31 will parameterise g + 1 independent cycles Cj = {(pi 6 [0,27r)} on the Jacobian 
along which the Poincare maps will be defined. Many of the techniques used in this section 
can be found in the book [3^] . 

2.1 Higher times and zero-curvature 

If one can rewrite the equations of motion of an integrable two-dimensional field theory in 
the form of a zero-curvature equation for a one-parameter family of 1-forms J(x), namely 

dJ(x) - J(x) A J(x) = 0, 

then this leads straight away to the construction of an infinite set of conserved charges by 
considering the parallel transporter Q(x) with connection J(x) along a loop winding once 
around the worldsheet. The flatness of J(x) immediately yields 

d a tr n(x) n = d T tr n(x) n = 0, VnG N. (2.1) 

Moreover, as was shown in [2] the invariants tr Q(x) n are in involution with respect to the 
Poisson bracket 

{tr n(x) n , tr n(x') m } = 0, Vn, m E N. (2.2) 

This condition contains (12. ip as a special case since the worldsheet energy S and momentum 
V are related to the leading order asymptotic of Q(x) near x — ±1. In fact ( \2.2h is the 
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statement of the invariance of tr Q(x') m with respect to an infinite family of higher flows 
generated by tr Q(x) n . We will now show that the Hamilton equations of motion corre- 
sponding to these higher conserved charges tiQ(x) n also take the form of a zero- curvature 
condition. 

Let us start by determining the evolution of the space component J\{x) of the lax connec- 
tion under the higher flows, namely {tr Q(x) n , Ji(x')}. For this we first obtain the following 
Poisson bracket with the transfer matrix T(<7i,<7 2 ,x) = Pexp JJ 1 dcrJi(a,x), 



{T((Ti, cr 2 , x)fJi(cr 3 , x')} = / da(T(a 1 ,a,x)®l){J 1 (a,x)^J 1 (a 3 ,x')}(T(a,a 2 ,x) 



(2.3) 

which requires the Poisson bracket {Jif Ji} given in first obtained by J.-M. Maillet in 
[22] in the context of the principal chiral model 

{Ji(«7,x)® Ji{a 3 ,x')} = [r(x,x'), Ji(a,x) ® 1 + 1 <g> Ji(a 3 , x')} 5(a - a 3 ) 

- [s(x,x'),J 1 ((T,x)®l-l®Ji(a 3 ,x')}5(a-a 3 ) (2.4) 

— 2s(x, x')5'(a — er 3 ), 

where 

.. 2vr x 2 + x' 2 - 2x 2 x' 2 ... 2tt x + x' . . 

rx, x ) = — 7 =- — — ^2-7/, s(z, x ) = — -=- — ^-7?. (2.5) 

VA {x — x')(l — ar)(l — x ) VA (1 — x ){\ — x' ) 

Inserting (12.41) into (12.31) . integrating by parts for the 5'-term and using identities like 

dT 

— (<Ji,a 2 ,x) = J 1 (a 1 ,x)T(a 1 ,a 2 ,x) 

^ (2.6) 
(cr 1 ,a 2 ,x) = -T(a 1 ,a 2 ,x)J 1 (a 2 ,x), 



da 2 
yields 

{T{<j u a 2 ,x)®Ji{<j 3 ,x')} 

= - 2(5(a 3 - a x ) - 5(a 3 - a 2 ))(T(a 1 ,a 3 ,x) g> l)s(x, x')(T(a 3 , a 2 , x) g> 1) 

+ e(ai - o- 2 )x(o- 3 ;(Ti,o- 2 )(T((Ti,(T3,x) <g> 1) (2.7) 
x [(r + s)(x,x'),J 1 (a 3 ,x) ® 1 + 1 <g> Ji(<r 3 , ^)K T (^3, 02, x) ® 1), 

where e(cr) = sign(cr) is the usual sign function and x(cr; 0"i, cr 2 ) is the characteristic function 
of the interval between <j\ and a 2 . If we are working on the circle, let 0\ = a + 2tt, a 2 = 
a, cr 3 = a' and identify the monodromy matrix as Q(a, x) = T(a + 2n, a, x) then the previous 
equation reduces to 

= (T(a + 2tt, a', x) ® l)[(r + s)(x, x'), Ji(a, x) ® 1 + 1 ® Ji(0, ^ x ) ® 
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Making use of f 1 2 . 6 j) again, we can rewrite this as 



{n(a,x)^J 1 (a , ,x')} = d a/ J(a,a',x,x') + [J(a,o',x,x , ),l® Ji(cr',x')} 



(2.8) 



where 



J(a, a',x,x') = (T(cr + 2ir,a' 1 x) ® l)(r + s)(x, x')(T(a', a, x) <g> 1). 
Taking the trace over the first factor of the tensor product yields 

{tiQ(a,x), Ji((j\x')} = d a ,J(a,a',x,x r ) + [J(a, a' , x, x'), ^(a' , x')], (2.9) 

where J {a, a', x, x') = tr x [(T(a', cr, x)T(a + 2n, a', x) <g> 1) (r + s) (x, x')] . In fact, using the 
translation invariance of the transfer matrix T by 2ir and the definition of Q(x), we see that 
J^(a, a', x, x') does not explicitly depend on a and can be written more succinctly as 



If we interpret the Poisson bracket {tr Q(a, x), Ji(<j', x')} in (12. 9p as the "time" derivative 
of Ji(a',x') with respect to the time generated by the Hamiltonian trf2(x) then (I2.9P takes 
exactly the form of a zero-curvature equation. This indicates that (I2.10p ought to be related 
to the Lax matrices corresponding to all the higher order flows generated by the Hamiltonians 
trf2(x), just as Jo and J\ were the Lax matrices generating r and a respectively. In fact, as 
we will show, one should Taylor expand (12. 9p and (12 . 1 0[) around x = ±1 thereby obtaining 
a discrete set of independent times t nj ±. 

An important remark is in order at this stage: since we are really doing string theory 
in conformal static gauge by imposing the Virasoro constraints and static gauge fixing con- 
ditions, which constitute a set of second class constraints in the Hamiltonian formalism, 
one should take care in imposing them consistently. This means that we should define an 
appropriate Dirac Bracket corresponding to every Poisson bracket and write everything in 
terms of those. Once this is done, the Virasoro constraints and static gauge fixing conditions 
can then be imposed without worry at any level of the calculation. However, as we show in 
appendix E2 for all the brackets of interest in the following, the Dirac and Poisson brackets 
are identical. Thus in the remainder of this section we shall denote brackets by {•, •} without 
specifying whether they are Dirac brackets or Poisson brackets. 

Let us first obtain the equations of motion for the monodromy matrix with respect to 
the Hamiltonian trf2(x). Starting from the Poisson algebra of the monodromies [221 E], 



J (a', x, x) = tn [(ft(</, x) <g> l)(r + s)(x, x')] . 



(2.10) 



{Q(x) f tt(x')} =[r(x, x'), Q(x) ® Q(x')} 

+ (Q(x) <g> 1) s{x, x) (1 ® n{x')) 
- (1 <g> tt(x')) s(x, x') (Q(x) ® 1) , 



(2.11) 



and taking the trace over the first factor of the tensor product as above yields 



{trfi(x),fi(x')} = [J{x,x'),Q(x')]. 



(2.12) 
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Once again, if we interpret the Poisson bracket {tr Q(x), Q(x')} as a time derivative, this 
last equation starts to take the form of the (a, r)-evolution equations 

[d T - J (x'), Q(x')] = 0, [d a - Ji(x'), n(x')] = 0. (2.13) 

The expression ( I2.10p for the Lax matrices can be simplified further. Using (12.51) the sum 
of the (r, s) -matrices entering in (I2.10p is 

/ \ / /\ 27T 2x 

r + s ){x,x) = ? =- — -77. 

VA [x - x')(l - x 2 ) 

Now by definition, t] = —t a <3t a where the su(2) generator t a is related to the Pauli matrices 
as t a = -^cr a . Thus ( 12.101) can be written as 

7r 2x 2 

J(a', x, x) = -=- —- — tr p(a', x) a a ] a a . (2.14) 

VA {x - x')(l - x 2 ) 

Now it is straightforward to show that for any matrix A G SL(2, C) the following is true 

v _ l t^jAaa\aa v = ^ where v _ lAV = ^ ^ A ^ 

— A_ 

i. e. V is the matrix of eigenvectors of A and A± are the eigenvalues. Since the eigenvalues 
of il(a',x) are e ±ip ^ and its matrix of eigenvectors is ^(x), this identity implies that the 
Lax matrix (12.141) corresponding to the Hamiltonian trQ(x) can be simplified as 

. 47ri sinn(x) ^/(x)a a ^(x)~ 1 

trflx — j( Xta f)= p y . 2.15 

^\ 1 — 1/ x z x — x' 

But since trf2(x) = 2cosp(x), it follows that the Lax matrix responsible for the flow of the 
Hamiltonian p(x) is 

p(x) <— > J(x,x) = -- 7 =^- T — — j . (2.16) 

Now we expand this around x = ±1 by extracting the Lax matrices associated with the 
Taylor coefficients of the quasi-momentum about x = ±1, namely 

res x=± i (x =F l)~ n p{x) < — > J n ,±{x) = res x=± i (x =F l)~ n J(x, x). (2.17) 

Using the straightforward identity for a rational matrix M(x) with singularities at x = ±1 

res x=±1 ^M = -(M(x')) ±1 , (2.18) 
x — x 
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where (M(x')) ±1 denotes the pole part of M(x') at x 1 — ±1, one can recast the Lax matrix 
( I2.17P in the much more useful form 

At this point we can also define the corresponding hierarchy of times t nj ± as the times 
generated by the Hamiltonians res x= -|-i (x =F l)~ n p(x) in (12.171) . namely we define 

d in ± = {res x=±1 (x T l)">(x), •} . (2.20) 

Going back to equation (12.121) . if we follow the prescription we just established to go from 
(I2.15P to (12.171) . namely of dividing through by —2 sinp(x) and taking the residue at x = ±1 
one readily finds the equation governing the evolution of the monodromy matrix under the 
hierarchy of times (12.201) 

[d in± -J n , ± (x'),n(x')] = 0, (2.21) 

which is exactly of the form (I2.13p . As an application of equation (I2.19P for the hierarchy 
of Lax matrices we show that the first two of these matrices Jo i= t are related to the original 
Lax connection J± = Jq ± J\. Indeed, applying the asymptotics for the quasi- momentum 
v( x ) ~x->±i — frf to equation (12.171) with n = we find 

~ ixi v&(±l)a 3 ^(±l)- 1 

" Jo ' ±(x) = ± 7x *Wl ' 

Now the components J± of the Lax connection are associated to a ± = |(t±ct) translations 



which are in turn generated by £ ± V = ^^ 2 ± and hence J±(x') = ^^Jo,±(x') since [T3 



2 ' 

vAre± J 



£±V <— > J 0)± (x') = -^ 7 M/(±l)a 3 ^(±l)- 1 = J ± (x'). (2.22) 

7T 1 =P X 

Finally we derive the evolution equations for the Lax matrices (12.191) under the hierarchy 
of times (12.201) and show that they take the zero-curvature form. We follow an argument given 
in [M] for finite-dimensional systems which applies readily here. Writing the monodromy 
matrix as Vl(x') = ^(x') di&g(e ip ( x \e~ tp ^) 1 &(x')~ 1 , equation (I2.2ip implies that 

^(x')" 1 (d in ± ^(x')) - ^(x')" 1 J n ,±{x')V{x'), diag(e ip(x) , e" ip(:c) )j = 0. (2.23) 

But any 2x2 matrix commuting with a diagonal matrix must itself be diagonal, and therefore 
we may write 

d in± V(x') = J n ,±(x')*(x') + ^(x')D(x'), (2.24) 
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for some unknown diagonal 2x2 matrix D(x'). Let us denote the multi-indices labelling the 
hierarchy, such as (n, +), using capital letters, e.g. N = (n, s n ) where n e N and s n = ±1. 
So let N = (n, s n ) and M = (m, s m ), then we have for Jn(x') = J n>Sn (x') 



dt M JN(x') 



Jm{x') 



2ni x' 2 ^(x')a^(x') 

in) 



(2.25) 



y/X x' 2 - 1 (x' 

where we have made use of (12.241) and the subscript on the commutator means we take the 
pole part of the whole commutator at x' = s n . Let us start by assuming that s n ^ s m , then 
Jm(x') is regular at x' = s n and only the pole part at x' = s n of the second term in the 
commutator contributes which is just Jjv(x'), so 

d t M JN(x') = [J M (x'),J N (x')] Sn , 

and likewise we also have &i N Jm(x') = [Jn(x'), JM(x')] Sm . Since [Jm(x'), Jn(x')] is rational 
with poles only at x' — ±1 and vanishes at x' = oo it can be written as a sum over its pole 
parts, namely 

[J M (x'), Jn(x')] = [J M (x'), Jn(x')}+i + [J M (x'), Jn(x')U. 

But because s n ^ s m we have {s m ,s n } = {±1} and the zero-curvature condition (12.261) 
below follows. If instead we assume that s n = s m , then we have 



~ 2m x' 2 ^(x')a^(x'y 1 j < , N 2ttz x' 2 ^(x / )a 3 ^(x')^ 1 
j n{x ) t=— t. — j-. — , Jm\x , 



\ in 







V\x' 2 -l {x>-s n y ' JUV ' J\x' 2 -l (x'-s n ) 

since both arguments in the commutator are regular at x' = s n = s m . The zero-curvature 
equation again readily follows from the above equation and (12.251) . i.e. 

dtJ N {x') - d~ tN J M {x') = [J M (x% Jn(x')). (2.26) 

Let us give an alternative basis J Uj ± for the string hierarchy whose zeroth level n = 
corresponds exactly to the Lax connection J±. If we define —7iK nt ± = res x= ±i(x =F l)~ n p(x) 
so that Ko,± = K ± then we have the following correspondence between integral of motion and 
Lax connection 

Vx _ VXf j , j \ 

In particular, from (12.221) we see that the zeroth level n = of this hierarchy is precisely 
the Lax connection J± associated with £ ± V = desired J ± = J±. It is 

straightforward to see by the linearity of the above expression for J n ^± and the constancy 
of the integrals of motion n n ,± that the new hierarchy is also commuting, namely it also 
satisfies the zero-curvature equation (12. 26ft . with d tn ± = {^K±K n ,±, •} 

d tM J N {x) - d tN J M (x') = [Jm(x'), Jn{x')}. (2.27a) 
Likewise, equation (I2.2ip also goes through unaltered and reads 

[d tM -J M (x'),n(x')] = 0. (2.27b) 
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2.2 Baker- Akhiezer vector and linearization 



Equations (I2.27P express the fact that the operators d tM — Jm{%') all commute among them- 
selves as well as individually with the monodromy matrix Q(x'). This means they can all 
be simultaneously diagonalised and there exists a solution if)(P') to the following equations, 
where P' = (x 1 , y') G V and Y : det (fl(x') — y'l) = is the spectral curve, 

[d tM - J M (x'))i,(P') = 0, VM 

(n(x')-y')^(P') = 0. 1 • } 

In this section it will be important to keep track of the explicit dependence of various 
functions on the hierarchy of times and so we will use the notation {t} for the complete set 
of times to,±>*i,=b • • • an d write for instance Jm(x', {£}), {£}) and i/j(P', {t}). 

The idea of finite-gap integration (see [131 E] and references therein) is to identify the 
analytic properties of the vector i/)(P', {t}) which specify it uniquely. To this aim we follow 
[TUGS] an d introduce the normalised eigenvector h(P', {£}) of Q(x') which is normalised by 
the condition a • h — 1 where ol = (1,1). Using this vector we can look for solutions to 
fl2T28|) in the form 

^(P>,{t}) = $(x',{t})h(P>,{0}), (2.29) 

where ty(x',{t}) is a formal matrix solution to (d tM — Jm(x'))^(x') = 0, VM so that 
(I2.29P trivially satisfies (d tM — Jm{x'))i^{P') = 0, VM. If we fix the initial condition to be 
i/?(P', {0}) = h(P', {0}) so that ^f(x', {0}) = 1 then by uniqueness of the solution with initial 
condition {0}) = tt(x', {0}) it follows that $(x', {t}){l(x', {0}) = n(x', {t})^(x', {t}) 
and therefore (12.291) is indeed also an eigenvector of Q(x', {t}). 

We now analyse the analytic properties of the vector i/j(P',{t}) in the form (I2.29P by 
obtaining the analytic properties of ^/(x') and h(P', {0}). First let us rewrite the hierarchy 
of Lax matrices in the more transparent form 

J n , ± (x') = (^(x')s n „ ± (x>3^(x , )^ 1 ) ±1 , 

where the singular parts s ni ±(x') are defined as 

v ± M= (-^(^+(^)) ±i - < 2 - 3 °) 

In the particular case of the zeroth level Lax matrix Jo t ±(x') the singular parts are precisely 
those of the Lax connection J± as defined in [T3J [Hj, namely so t ±(x') = -^j- Because 
J n ,±{x') only has poles at x' = ±1 it follows by Poincare's theorem on holomorphic differential 
equations that ty(x') is holomorphic outside x' = ±1. By studying the asymptotics of the 
equation for ^(x 1 ), its behaviour near x' — ±1 is easily show to be 

$(x', {t})e~ E « s "-±'"- ±fT3 = 0(1) as x -> ±1, 
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where 0(1) denotes a matrix holomorphic in a neighbourhood of x' — ±1. Moreover, using 
the fact that </„, i= t(oo) = we observe that <9f M \l/(oo, {t}) = 0, VM and hence \l/(oo, {t}) = 1 
by the choice of initial conditions. Turning to the normalised eigenvector h(P', {£}), a 
standard analysis of its analytic behaviour reveals that it is meromorphic in P' and uniquely 
specified by the following condition 

(h) > 7(W) _1 oo-, h 1 (po + ) = l, and (h 2 ) > 7({t})~ 1 oo + , ^(ocT) = 1, 

where the divisor 7({t}) of degree g + 1 is called the dynamical divisor. The analytic data 
gathered above for ^(x') and h(P', {0}) is sufficient to uniquely characterise the components 
of i^(P ; , {t}) as Baker-Akhiezer functions, namely 

C0i) > 7cT locr 7 ^i(oo + ) = 1, and (ip 2 ) > 7 c T 1 oo + , ^2(00") = 1, 

UiOA {t})e TE " s "- + '"'+ = 0(1), as x -> 1, 
with < „ 

{^i{x ,± , {t})e TE "^-*™- = 0(1), as x -1, 

where 70 = 7({0}) is the initial divisor. Notice that the hierarchy of times enters linearly 
in the definition of the Baker-Akhiezer vector i/)(P', {t}) through the essential singularity. 
This is a very general feature of finite-gap integration. When explicitly reconstructing the 
Baker-Akhiezer vector satisfying the above conditions in terms of Riemann ^-functions on 
S, the singular parts give rise to a unique normalised Abelian differential of the second kind 
dQ with poles at x' = ±1 of the prescribed form 

r ir\ -jo / 11 1 f S + (x'^, {t}) = ± ^2 s n.+ { x ')t n ,+ , 

dQ = idS±, asx ±1, where <^ /± ) \( ^ n ^} ' ^ 

\ S^{x ± ,{t}) = ±}^ n s n ^(x)t n: -. 

All the time dependence of the Baker-Akhiezer vector, and hence of the solution, is encoded 
in this meromorphic differential dQ which is linear in the hierarchy of times. In fact, we can 
define a differential associated to each time of the hierarchy by writing 

dQ = t n,+d^ n ,+ + Yl t^-d^n,- = ^ t N dQ N , (2.31) 

n n N 

using the multi-index notation, where the normalised Abelian differentials of the second kind 
dQ n ,± are defined uniquely by their respective behaviours at the points x' = ±1, namely 

d^l ni+ (x' ± ) = ±ids rii+ (x') as x' — > +1, d^l n ^(x' ± ) = ±ids n _(x') as x' — > —1. 

This correspondence between times of the hierarchy and meromorphic differentials on E 

t n ,± ^ d£l nt ± 

is a very general feature of finite-gap integration. We say that the differential couples to the 
time for obvious reasons from (12.311) . As we saw in the previous sections, every Hamlitonian 
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corresponds to a Lax matrix which is responsible for generating the corresponding time in 
the Lax formalism. Here we see that every Hamiltonian also corresponds to a meromorphic 
differential on E responsible for generating the corresponding time in the finite-gap language. 
Notice the splitting between differentials singular at x' = +1 and those singular at x' — — 1. 
These are related to left and right movers of the string. For instance, at the zeroth level 
n = we have a ± = = — i ,± and dq± = dq ± dp = —2iidQ 0j ±, so in particular 

t 0+ dQ, + + t dtt = -^-(crdp + rdq), 

which is the usual dQ defined in [13l HU [15] where all the higher times are set to zero. 

In the next section we will be perturbing finite-gap solutions and so we give here the 
explicit formulae for the generic finite-gap solution in terms of Riemann ^-functions on E. 
Details can be found in [131 El D3] • Of particular interest for constructing the embedding g 
of the string in SU (2) 



* = { -i ii ) e su ^- 

is the dual Baker-Akhiezer vector which is defined relative to the conjugate divisor f7o and 
opposite singular parts —s nt ±. The components are explicitly constructed as [TH [T5] 



Z x = Cy>+(0 + ), Z 2 = C , ^ 2 + (0 + ), (2.32a) 

where C G 1R is a normalisation constant chosen such that \Zi\ 2 + \Z 2 \ 2 = 1 and x{P) is a 
meromorphic function on E with divisor (x) = 7 • t^q ■ B~ l and normalised by x{°° + ) = 1 
(B is the divisor of branch points of E). The components of the dual Baker-Akhiezer vector 
at + are explicitly given by 

^+0+ = h_ 0+ L-LJ-V — J°o± Jb^L ^exp +- / dQ 

^ 0+ = h + + )/ ^ Jo ° w % ^-exp -- / dQ 

9(f b dQ + D;U)e(27r £_ W -£> ; n) \ 2 





2.3 Quasi-actions 

Remember that the Lax matrix in (12.161) is responsible for the flow of the Hamiltonian 
trQ(x) = 2cosp(x). Thus going back to the corresponding Hamilton equation in Lax form 
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we can rewrite it as 
y/X 



2m 



]7i 2 i 



1 - 



x- 



p(x), Ji(x') 



d a - Ji (x 



/N ^/(x)^(x)- 1 



x — X 



(2.33) 



Integrating this equation in x over the different a-cycles, and recalling that the action vari- 
ables are defined as Si = J a (l — ^2) p(x)dx we find 



{Si, Ji(x')} 



d a - JAx') 



Air 



*(x)cr 3 *(x) 



-dx 



x — x 



(2.34a) 



and similarly integrating around the point x = 00 and recalling that the global SU(2)r 
charge is defined as ^ = ^7 § (l — -^)p{x)dx we find 



\{R,JiW)} 



-dx 



x — x 



(2.34b) 



Equations (12.341) simply say that the Hamiltonian flow of the action variables Si and R are 
generated by the following respective Lax matrices 



Si 

R 
~2 



M x ) = t ; — dx > 

47T /„ . X — X' 



An 



^/(x)a 3 ^/(x)~ 1 



(2.35) 



x — x 1 



dx. 



Because any integral of motion can be expressed in terms of the action variables Si, one 
ought to be able to use equation ( 12. 35ft to derive the Lax matrix for any other integral of 
motion. Indeed, for instance we know that 



i=i ^ 
9 



8=1 

and so this means one can write 

9 



{S±V,-} = J2 



i=l 



2n 

dq 
2n~ 



dq± 
2tt 



^ Isr. 



00 
00" 



2tt/ 2 
1 



(2.36) 



2tx 2 



5R 



{Si,-} 



dq±\ 1 



2tt J 2 



{R,-} 



where dq± = dq ± dp. Making use of the Lax matrix for the action variables (12. 35ft and the 
fact that the differentials dq± are normalised as j a dq± = means we can write the Lax 
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matrix for S ± V as follows 



1 9 

47T ^ 



^(x)^*^)- 1 ^ f dq±_ f ^(x)a3^(xy 1 dx f dq± 



X — ./&• 27T y^,. X — X J n . 27T 

1 i ^(x)a^{x)-\ dx f°° + dq ± 



i=l 

An x — x' ./™- 27r 

Written in this form we can apply the Riemann bilinear identity to obtain 

£±P < — ► -2(res I= i +res x= _i j — — dx, (2.37) 

V /XX 

where an overall factor of two came from the fact that we get equivalent contributions from 
both sheets, namely at x ± = (+1) ± and x = (— 1) ± . Note also importantly that there is 
no contribution from the apparent pole at x — x'. This is because x = x' is not actually a 
pole of the Lax equation itself, as can be seen from f!2.33j) which is perfectly regular as x 
approaches x' since 

trivial fact that diagonal matrices commute. An equation such as (I2.37P relating an integral 
of motion to a Lax matrix should really always be understood as a relation between two 
ingredients of a Lax equation. To evaluate the residues in ( 12.371) we note that the Abelian 
integrals q±(x) have poles at x = ±1 with the following asymptotics 

2nn + 27tk_ 



q+{x) ~ x _ +1 -, q-{x) 



x — 1 X + 1 

It follows now using the identity (I2.18P that 

e ± v — (^(xyaMxr 1 ^-) = ^ 7 *(±i)a 3 *(±i)- 1 = j ± (x'), 

and we recover exactly the same expression as before (12.221) . It is important to note that it 
was the multivaluedness of the Abelian integral q± (P) = J P dq± (or equivalently the fact that 
dq± had some non-trivial periods) which resulted in a non-zero answer for the corresponding 
Lax matrix. Indeed, the Lax matrix obtained by this argument clearly depends only on the 
cohomology class [dq±] G i? 1 (S, oo^) of the Abelian differential dq± one starts off with on the 
singular algebraic curve H/{oo ± }. One can see this explicitly from the equation preceding 
(I2.37P or otherwise from (12.371) itself: suppose dq±,dq' ± are two representatives of the same 
cohomology class, then dq± — dq'± = df is exact and the corresponding difference of the 
expressions in (12.37D is 

j_ T vp)*^)- 1 m dx 

2 4^ + x(P) -x' 2vr 

p=(±i)± v ' 
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where ^(P) = (h(P),h(aP)). But this is the sum over the residues of a well defined 
meromorphic differential on T,/{oo ± } (since f(P) is single-valued) and so is zero. 

One could use the same trick as above to compute more explicitly the Lax matrix for the 
action variables ( 12.35(1 . To simplify the notation we first combine Si and \R into the g + 1 
filling fractions 

V* f A 1 



Sl = -^Hj Al [ 1 -^) p{x)dx 

where Ai is the cycle going around the I th cut counterclockwise on the top sheet. They 
satisfy JTJ/ii = \{L — R) where L is the global SU(2) L charge. So to apply the previous 
reasoning we could write 

3+1 

j=i 

For the same argument to follow through we must introduce Abelian differentials dq^ of 
the second kind (so dq^ has no residues) such that 

dq^ = 0, [ dq^ = 5 u , (2.38) 
'Ai JBi 

where Bj is the contour going from oo + to oo~ through the 7 th cut. Such differentials exist: 
consider g + 1 independent differentials from the hierarchy, and call them dQj. Then the 
(g + 1) x (g + 1) matrix Ajj — f B dVLj is invertible, and so dq^ = A~^jdVL K have the desired 
property. Yet since the conditions (12.381) on the differentials dq^ uniquely specify their 
cohomology class in i? 1 (S, oo^), by the preceding remark they are also sufficient to uniquely 
fix the resulting Lax matrix. By the procedure of section [2721 these Lax matrices yield unique 
normalised Abelian differentials which satisfy ( 12.38(1 . which we still denote dq^ by abuse of 
notation. Since the operations of constructing a Lax matrix from a given integral of motion 
and that of constructing an Abelian differential from a given Lax matrix are both linear, it 
follows that the equation for 8 in (I2.36j) translates into an equation in terms of differential 
forms on ^/{oo^}. Rewrite this equation as 

- - 1 ( 1 1 - r + ( r *[ k*+±*) - e ( l 

(2.39) 



i=i 



9+1 



it follows that 

d 1=X,{J B d l) d 1 (1) i ( 2 - 4 °) 
and in particular this leads to the following equation which will be important later 

9+1 



/ <*9 = £(7 d( i) I d * {1) 

J Bq I=l KJbj J Jbq 
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where Q G T is a singular point on T which is blown up on the desingularised curve E to 
two points Q G X and Bq is a curve joining Q having zero intersection number with any 
of the a- or 6-cycles, see Figure ITD 




Figure 11: Definition of the cycle Bq for a given singular point Q. 



3 Perturbations of finite-gap strings 

Given a finite-gap solution Z\ with underlying algebraic curve £ of genus g, one can obtain 
its stability angles by considering nearby solutions Zi + 5Zi with algebraic curves S £ of 
genus (7 + I. In other words, perturbations of a given finite-gap solution Zi correspond to 
degenerations of a genus g + 1 algebraic curve S € into the genus g curve £ of the solution 
Zj, see Figure [T2J Now since we are concerned with real finite-gap solutions, constructed 




Figure 12: Pinching an a-cycle. 



from real algebraic curves E (see [T3] for a discussion of reality conditions), the degeneration 
process in Figure [12] describing the perturbation should respect this reality condition. This 
forces us to consider degenerations through the pinching of imaginary cycles, which we can 
choose to call the a-cycles as in [13] . We discuss the pinching of a-cycles in appendix [D] 

As discussed in section 12.21 the dependence of the general finite-gap solution on the 
hierarchy of times {t} is entirely encoded in the normalised Abelian differential of the second 
kind dQ = ^2 N t]ydQN defined in (12.311) which enters the reconstruction formula as follows 

z ' = Ci 6U b dQ + D] n) exp J Pi da ) • (31) 

where P± = oo + and P2 = oo~. In this expression we have hidden all the time independent 
part into the overall constants Ci whose specific forms can be retrieved from the complete 
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reconstruction formulae (12.321) . A nearby solution Z± + 5Zi is constructed with the same 
formulae but from slightly deformed data (which includes a deformed curve E e ) 



Zi + 5Zi = CI 




exp — i 




o+ 



) 



(3.2) 



The ingredients of the deformed solution are as follows. First of all, since the underlying 
curve E e has genus g+1, the arguments of the ^-functions for this curve are (g+l)-component 
vectors, namely 



In the singular limit e — » one has b e — > b and uj 6 — > u which are the 6-cycles and the g 
holomorphic differentials on E respectively. The extra 6-cycle b e becomes a degenerate cycle 
on the curve E, see Figure [TJJ In appendix Owe show that the extra holomorphic differential 
cUq on E e acquires a simple pole at the singular point and so becomes a normalised Abelian 
differential of the third kind. The Abelian differential dQ e on E e is defined by the same 
singular parts (I2.30P as dQ at x — ±1 but could potentially acquire an extra simple pole at 
the singular point. However, because dQ e is normalised on E e , its residue there would vanish 
in the e — > limit, so that in fact dQ e — > dQ. One can also show that C\ — > Ci as e — > 0. 

The important object in (13.21) when considering the singular limit e — > is the period 
matrix which can be broken down into blocks in a natural way 



The singular limits of each block follow from the above considerations of b e ,u3 e in the limit 
(see appendix |D1 for details). In particular, IT — > IT as e — > which is simply the period 
matrix of E. The vectors IIq also stay finite in the limit. The top left component IIq on 
the other hand diverges in this limit, leading to a simplification of the Riemann ^-function 
$(■; IT) as e — >• which becomes expressible in terms of the Riemann ^-function #(•; II) of E 



Now taking into account all the above limits and working to first order in e, a direct 
computation shows that the difference 8Z^ between expressions ( 13. 2 j) and (13.11) contains 
three types of contribution 



5Zi = ( {periodic} + {periodic} x e ifb o dQ + {periodic} x e' ifb o dQ ) x e" n o°, (3.3) 



where "{periodic}" denotes functions periodic in all the angle variables tpj of the underlying 
finite-gap solution ( 13. lft , namely invariant under tpx — > <Pi+ 27r for each 1 = 1,.. . ,g + l. The 



D 





as in (jH5j) . 



34 



three contributions in (I3.3P correspond to three different stability angles of the underlying 
solution (13.11) which can be read off directly 



0, v { ±> =±2tt dq (1) , I = l,...,g + 1. (3.4) 

</f>o 

The zero stability angles are related to the (/^-translation invariance of the equations of 
motion which is explicitly broken by the finite-gap solution (13. ip . 

Now stability angles are defined modulo 2n but for the underlying solution (13.11) to be 
periodic requires that 

pOO + 

2tt / dq {1) E 2ttZ, / = l,...,g + l, 

J oo~ 

therefore we can redefine the stability angles as 

< (/) - ±2tt ( / dqW + r dqW) = ±2tt / dq™, (3.5) 

\Jb J 00+ I J B 



where the contour Bo runs from oo + on the top sheet to oo~ on the bottom sheet, by going 
through the th cut, see Figure [T3"l In the singular limit e — > the th cut shrinks to a point, 
say P and so (13. 5p yields 



(ffo poo \ I pfo poo 

\ dq {1) + \ dq {l) = ±2tt / dq {1) - / a*dq {1) 
Joo+ J Pq J \Joo+ J Pq 

(fPo poo + \ / rPo pa-Pa \ 

/ dqW - / dg« = ± 27 r ( / dqW + / dq^ ) 
Joo+ J aPo J \Joo+ Joo+ J 



(3.6) 



'vP 

= ±2tt ( g «(P ) + q {I \aP )) = ±Anq^(P ), 

where q^\P) = dq^ with the integral running along the top sheet (the precise choice 
of contour then doesn't matter since dq^ is normalised) and the last equality follows from 
P = &P by virtue of P being a singular point. 

By repeating the calculation in (13. 6p but for the £> -period of dp (the integrality of the 
fe-periods of dp follows from the closed string requirement, namely that the finite-gap solution 
be periodic under a — > a + 2%), 

/ dp = 27m , no G Z, 

•/Bo 

the details of which can be found in [T3], we arrive at an equation for the location of the 
singular point Po, namely 

p(P ) = n n. 
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Figure 13: The canonical cycles before (a) and after (b) shrinking of the th cut. Note that 
it doesn't matter where the shrinking cut lies with respect to the other cuts, but for the sake 
of clarity of the figure we chose it to be the furthest to the left. 

The above analysis shows that to this singular point Po there corresponds two stability angles 
for each of the g + 1 cuts determined by the So-period of corresponding quasi-action dq^ 1 ' or 

4 ] = ±4W / )(P ). (3.7) 



4 Semi-classical energy spectrum 



As we recalled in section 12.21 every finite-gap solution to the equations of motion of a 
bosonic string on M x S 3 is constructed from a finite-genus algebraic curve S equipped with 
an additional set 70 of g + 1 points on it called a divisor (of degree deg % — g + 1). This 
algebro-geometric data can be identified with a bundle Ai^ g+2 ^ over the moduli space C of 
the algebraic curve S, of dimension dime C = g + 1, 

S 9+1 (Z) -> M<£ 9+2) -> c, 

whose fibre over every point of the base, specifying a curve S, is the (g + l)-st symmetric 
product S 9+1 (E) = Hi 9+1 /S g+ i of £ (see [13] for moer details). The finite-gap construction 
of [T3] defines an injective geometric map [HH] from this algebro-geometric data Ai^ 9+2 ^ into 
the space of complexified solutions j e sl(2, C) to the equations of motion of a string 
moving on R x S 3 which also satisfy the Virasoro and static gauge conditions, 

Q : A4 29+2) - Si. (4.1) 

Since a general point in phase-space is the restriction to the hypersurface r = of the general 
solution we can identify the space Sc of (complexified) solutions with (complexified) phase- 
space . Furthermore, the subset Sq C Sc of solutions satisfying Virasoro and static 
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gauge conditions can be identified with the second class constraint surface Vq C defined 
by these conditions. We can describe the map (14. 1H as an embedding M (2g+2) ^V^. If we 
further impose reality conditions by restricting the algebro-geometric data Ai^ 9 ^ to real 
algebro-geometric data (see [13] for a detailed discussion of reality conditions) then finite-gap 
integration describes an injective map [35] 

& : >4 29+2) ^ PL (4-2) 

from the (g + l)-dimensional toric fibration T 9+1 — > A4^ 9+2 ^ — > £r, with dim^ £r = g + 1, 
into the (real) phase-space of strings on txS 3 satisfying the Virasoro and static gauge 
constraints. Introducing the inclusion by : <—* of the second class constraint surface 
C Vr, the Dirac bracket on is the pull-back of the symplectic structure u> on V^. 
As was show in [2], the pull-back to M.^ 9 ^ of this symplectic structure l v uj on takes 
the simple form 

9+1 
1=1 

The different variables in (14.31) are defined as follows 



The action variables Si are given by the filling fractions 



Si = tt^t / a, I = l,...,g+1, 



2m 



where Ai is the cycle encircling the 7 th cut Ci on the physical sheet of £ represented 
as a hyperelliptic curve and a = ^zalp is a special 1-form on £, with z = x + \ and 
p(x) being the quasi-momentum. 

The angle variables tpi are specified by the image of the divisor 70 on the generalised 
Jacobian J(£, oc^) of the curve £ under the extended Abel map A : S 9+1 (T,) — > 
</(£, oo 1 * 1 ), or more precisely 

<Pi = Ai(%) - A g+ i(%), i = l,...,g, ^9+1 = ~A+i(7o)- 



The injective map ( 14.2ft can thus be thought of as an embedding in phase-space of a (<? + 1)- 
parameter family of isotropic ((7 + l)-torii parametrised by {SV}?=i since the pull-back ( 14.31) of 
the symplectic form uj to these torii is identically zero. This is the necessary set-up to apply 
the Bohr- Sommerf eld conditions (11.121) for the quantisation of a p-torus in an n-dimensional 
phase-space, where here the total phase-space is infinite dimensional so that n = 00 and 
p = g + 1. The condition ( 11.121) also involves the stability angles of perturbations around 
the p-torus which we computed in the section [3j So applying ( 11.121) to the finite-gap string 
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we can write down the Bohr-Sommerfeld quantisation conditions for the action variables of 
the string as follows 



f^+i+E („.+!) £ +0( »), (4.4, 

where we have used the fact that the Maslov index for the ^4/-cycle (J = 1, . . . , g + 1) in the 
generalised Jacobian J(S,oo ± ) is simply fij = 2. We emphasise that (14.41) is only valid in 
the harmonic oscillator approximation Ni ^> n a where the perturbations are much smaller 
than the background filling fractions. 



4.1 The main result 

In the semiclassical regime, the Hamiltonian is defined by the same classical function of 
the actions E c \[Si, . . . , S g+ i] but evaluated on the action operators since by semiclassical 
integrability we have that [Si, Sj] = 0(ft?), so 

"^string = E c \[Si, . . . , Sg+l] + 0(h 2 ). 

It follows that the energy spectrum is simply the classical energy E c \ evaluated on the 
eigenvalues of the action variables (14. 4p namely 



E = E, 



cl 



Nih + ^+ ( n <* + l) 



2 2tt 



a=g+2 



(9+1) 

n 

2V 



-h 



0(h 2 ). 



We now Taylor expand this using the fact that Nj ^> n a to obtain 



E = R 



cl 



iVi 



%,...,[ N g+l + 



h 



Using equations (12.391) and (13. 5p to express dE c \/dSi and v& respectively as £>-periods 



9+1 oo 

+ EE 

I=\ a=g+2 
(I) 



1 \ dE d V a 



n r 



(I) 



2 J dS! 2ir 



h. 



E = E, 



cl 



h,...A N g+l 



h 



g+l co 

E E 

i=l a=g+2 



2^ 



dg (/) /i. 



So 



where £> Q is the contour running from oo + to the singular point labelled a on the top sheet, 
and back on the bottom sheet to oo~. The sum over / can now be performed using equation 
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(g30D which yields 



E — E. 



cl 



N 1 + -)h,...A N q+1 + 



+ £ 



n a + 



2tt' 



/?. 



If we now formally think of the function E c \ as depending on the infinite set of filling fractions 
{<Sj}?=i) {^a}^=g + 2 but finitely many of which are turned off for the classical finite-gap 
solutions) then we can interpret the £>o,-period of dq/2ir as dE c \/dS a using a formal analogue 
of (12.391) for an infinite gap solution. One can then resum the resulting Taylor expansion to 
obtain the following formal expression for the semiclassical energy spectrum 



E = E r 



N 1 + -)H,..., 



N, 



3+1 



n g+ 2 + - I h, 



(4.5) 



We stress that this is only a formal derivation as rigorously one would have to regularise the 
divergent infinite sum over stability angles at the intermediate steps as well as subtract off 
the energy of the vacuum (i.e. the zero cut finite-gap solution). But formally at least the 
result of the above derivation is the following: 

• The semiclassical energy spectrum is obtained by evaluating the classical energy func- 
tion of an infinite-gap solution on filling fractions quantised to half-integer multiples 
of H. 

• The infinite number of singular points of the spectral curve det (Q(x) — yl) = which 
accumulate at x — ±1 must be filled with half a unit of h in their ground state with 
an additional integer multiple of h for excitations. 



4.2 Comparison with alternative approach 

In [16] an alternative method was proposed for extracting the semiclassical energy spacing 
around any given classical solution from the algebraic curve E itself, without making use of 
the divisor 70 on E, and which the subsequent papers [TTJ, [18] built upon. The heart of the 
method resides in the assumption that the filling fractions Sj become quantised in integer 
units at least in a semiclassical approximation. This assumption seems natural because 
the filling fractions constitute the action variables of the theory (a fact proved only in the 
1R x S 3 subsector [H]) and we expect that after Bohr-Sommerfeld quantisation the action 
variables become half-integer multiples of h. Yet we see from (14.41) that this is not the 
case. As we have argued at the start of this section, a finite-gap solution can be pictured 
as a degenerate isotropic (g + l)-torus within the full infinite dimensional phase-space. And 
although the Bohr-Sommerfeld conditions (11.91) for an integrable system do imply that the 
action variables become half-integer multiples of h, we saw in section 11.21 and appendix 
IB1 that these conditions receive 1-loop corrections, when applied to a degenerate isotropic 
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torus, from fluctuations transverse to the torus in the form of stability angles. Only after the 
calculation in the previous subsection can one conclude that the semiclassical spectrum of 
a finite-gap solution is given by the classical energy of an infinite-gap solution whose fillings 
are half-integrer multiples of h. 

Now the algebraic curve £ is characterised by the quasi-momentum p[x) used to define 
the filling fractions Si as 

S I = 7 L^[ ( x + -)dp=-^-^ [ (l-^-]p{x)dx, 7=1,. ..,0+1. (4.6) 



2ni An J Ai \ x J 2ni An J A \ x 2 r 

The integer quantisation of these filling fractions in the semiclassical limit can be interpreted 
in the language of the gauge theory side by attributing to a single Bethe root one unit of 
filling fraction. In the semiclassical quantisation of a solution each cut of its algebraic curve 
thus turns into a large clump of Bethe roots with the filling fraction counting the number 
of such roots [TB]. The idea of [TUJ for obtaining the semiclassical energy spacings is then to 
compare the energies of two neighbouring classical solutions differing only by a single Bethe 
root. If the underlying solution is characterised by the quasi-momentum p(x) and has K 
cuts Cj with mode numbers rij, j = 1, . . . , K, 

p(x + z'O) +p(x - z'O) = 2nfij, x 6 Cj, j = 1, . . . , K, (4.7) 

then its perturbation is characterised by a perturbed quasi-momentum p(x) + 5p(x) with 
still the same K cuts but also with an extra isolated Bethe root at xk+i with mode number 
n K +i 

p(x + zO) + 5p(x + iO) +p(x - iO) + Sp(x - iO) = 2nn j , x E Cj,j = l,...,K, (4.8a) 

p(x K +i) + Sp(x K+1 ) + p{x K+ i) + Sp(x K+1 ) = 2nn K+1 . (4.8b) 
By using (14.71) we may simplify (j4.8aj) to 

Sp(x + zO) + 5p(x - iO) = 0, x e C j: j = 1, . . .,K. (4.9a) 

and since 5p(x) is small, by working to lowest order we can approximate (14.8bl) as 

pipK+x) = im>K+i, (4.9b) 

Equations f)4.9p are the starting point in [TU] for obtaining the semiclassical energy spacings 
by reading them off from Sp(x). 

Let us now show that the semiclassical energy spacings obtained by this method agrees 
with the semiclassical spectrum ( 14. 5ft obtained in the previous subsection. We know from 
( I2.39j) that the variation of the energy E of a classical solution as we vary the moduli Si is 



<?+i 



i=i 



dq 

Bi 
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It follows that adding a single Bethe root (which would correspond to setting 5Sj = h for 
some J) should increase the energy of the solution by 




(4.10) 



This is exactly what one gets if we set Nj —>■ Nj + 1 in (14 .5p and Taylor expand in the 
J th entry using Nj ^> 1. We easily find that the energy evaluated on the solution with 
Sj = Nj + 1 is equal to the energy evaluated on the solution with Sj = Nj plus the 
perturbation (I4.10p . Thus (14. 5 p predicts the same energy spacing (14.101) as we would expect 
if Bethe roots carried h units of filling fraction. 

Note finally that the energy E c \ we have been using is not the space-time energy of the 
classical solution but rather the worldsheet energy or the Hamiltonian of the fields Zi in the 
action. It can however be related to the space-time energy A by the following simple formula 
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5 Summary and Outlook 

We have obtained the semiclassical energy spectrum of bosonic string theory on 1 x S 3 as 
expressed in equation (14. 5 p by semiclassically quantising the general finite-gap solution of 
this theory. The derivation of H4.5|) can be summarised as follows. We have argued that the 
generic finite-gap solution can be thought of as an embedding of a (g+l)-torus into the full 
infinite dimensional phase-space of the theory. Since these torii are finite-dimensional, they 
are all degenerate isotropic torii located on the boundary dS of the infinite region S = {Si > 
0, V/}. But a procedure due to Voros [221 ESj provides a way of semiclassically quantising 
such degenerate torii: the method consists of studying neighbouring orbits in the small 
oscillator approximation, which would live on a neighbouring non-degenerate torus in the 
interior of S, and then quantise this torus in the usual way using Bohr-Sommerfeld-Maslov 
quantisation conditions. The computation in section 14.11 consisted in formally rewriting the 
quantised energy of such a linearised torus in terms of the energy of the infinite-gap solution 
it is approximating in the interior of S (but still near dS). The result is that the semiclassical 
energy spectrum can be obtained by evaluating the classical energy function on points of the 
following infinite lattice in S, 

Srehf^ + nYw. (5.1) 

Yet because we computed the semiclassical spectrum around finite-gap solutions which only 
describe dS, it follows that (14. 5 p only describes this lattice structure near the boundary 
of S where the n's are much smaller than the finitely many iV's in (14.51) . However, since 
the number of iV's is finite but arbitrary, by formally considering the infinite genus limit of 
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finite-gap solutions, the complete spectrum is described by E = E d [(Ni + |) h, . . .] so that 
( 15.11) should give the correct lattice structure in the whole bulk of S. 

Such a procedure for semiclassically quantising finite-gap solutions should be sufficiently 
general to apply with little modification to more general settings and in particular to the 
case of superstrings on AdS$ x S 5 . As stated in the introduction, it would therefore be very 
interesting to obtain the divisor for the full algebraic curve of AdS^ x S 5 by constructing the 
finite-gap solution in full generality on AdS 5 x S 5 . 

Finally, in view of ultimately obtaining an exact quantisation of string theory on AdS§ x 
S 5 we have argued that operator ordering issues will be of crucial importance since they 
already appear in the semiclassical analysis. In this paper we assumed for simplicity that 
the cohomology class of the subprincipal form vanished since with this assumption we were 
able to reproduce the semiclassical spectrum of [HI [T7J [18] at least for the fluctuations in 
the R x S* 3 subspace. This rules out many operator orderings for the exact quantisation, 
namely all those for which the action variables have a subprincipal Weyl symbol. 
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A Symbolic calculus of Pseudo-differential operators 

The passage from a classical system on phase-space T*X to its quantum counterpart involves 
promoting the algebra of classical observable C(T*X) to a noncommutative algebra A of 
operators. Classically, the Poisson algebra of observables is uniquely specified by the choice 
of a symplectic structure u = £\ dxi A d^ and the Poisson bracket of two observables f,g£ 
C(T*X) is then defined by {/, g} = uj(Xf,X g ), where Xu denotes the Hamiltonian vector 
field associated to any function H G C{T*X) satisfying ix H u = dH. To pass to quantum 
mechanics, the prescription of canonical quantisation is to promote the special functions 
Xi, £j G C(T*X) to operators x iy £j and the symplectic structure ui = dxiAd^i to the Weyl 
algebra = ihSy which admits the unique representation £j = = —ikd/dxi = 

—ihdi in terms of differential operators on L 2 (X). The problem that remains after canonical 
quantisation is to associate with any other given observable / G C(T*X) (which is a function 
of Xi, £j) a (pseudo-) differential operator / on L 2 (X), and it is immediately obvious that this 
is by no means unique. Many different operators correspond to the same classical function: 
for instance, given any t G R, the differential operator txidi + (1 — t)d\ ■ xi is a possible 
candidate for the quantisation of the function xi£±. In other words, it is not possible to 
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specify the operator ordering in an operator / starting from just single function / G C(T*X). 
However, with an infinite set of functions fk G C{T*X) it turns out to be possible to associate 
a unique operator / by canonical quantisation. Such a set defines a function of h through 
the asymptotic expansion 

f(z,bn)~52fk(x,Z)ti e . (A.i) 

k>0 

We refer to such a /i-dependent function f(h) G C(T*X) as a classical (Weyl) symbol, which 
is technically required to satisfy certain estimates, such as all its partial derivatives being 
uniformly bounded by some order function. 

Without going into details of the construction, we now state the map from symbols 
to pseudo- differential operator^ (\l/DO for short). Given a symbol f(h), we define the 
corresponding ^DO by specifying its action on u G L 2 (X) using the Weyl quantisation 
formula [30] 

(Op]T(/(ft))«) (x) = J e *(-»)-e/ (^±y, £; h\ u(y)dyd£. 

It is important to note here that the choice of Weyl quantisation in the definition of the \l/DO 
from its symbol does not limit us to having only Weyl ordered \&DOs. Indeed, the operator 
Op^(f(h)) is Weyl ordered only when the corresponding Weyl symbol is fo- independent. So 
it is precisely the subleading terms in the asymptotic expansion fTA~T]) of the symbol f(x, f ; K) 
which account for the different possible choices of orderings in the definition of the \l/DO. 
For example, the Weyl ordered operator of the classical observable is given simply by 
the Weyl symbol Xi£\, namely 

—ih 

OpJT(zi6i) = -y- ( x i d i + d i ■ x i) > 

whereas the left ordered operator —ihx\d\ which corresponds to the same classical observable 
Xi£i as Op ft (xi£i) is given by a Weyl symbol with a subleading term in h since 

Op^ + yj = -ihxidi. 

Naturally the right ordered operator —ihdi -X\ has Weyl symbol X\^\ — ~. A general \l/DO A 
always has a unique Weyl symbol, which is a ^-dependent function f(x, £; h) denoted a w (A). 
The leading non-zero term in the asymptotic expansion 1 1A.1I) of this Weyl symbol is called the 



principal symbol, denoted a^(A), and the subleading term is called the subprincipal symbol, 
denoted a^ h (A). For instance, if f (x, g) ^ then a^(A) = f (x, f ) and a^ h (A) = f x (x, £)h. 



5 Whcn the symbol f(x,£;H) is a polynomial in x,(; the associated operator is an ordinary partial differ- 
ential operator. To include the more general case when /(x, £; h) might not be a polynomial we talk about 
pseudo-differential operators. 
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An important object for the study of quantum integrability is the commutator [A, B] of 
two operators A and B. In the present context of \I/DOs one can show that if A, B are \l/DOs 
then their commutator [A, B] is also a ^DO with principal symbol 



(so that —iha t 




B Bohr-Sommerfeld for isolated periodic orbit 

Let 7 be a given periodic orbit of energy E, i.e. 7 C We henceforth assume that E is 
a regular value of H so that is a smooth codimension one submanifold of T*X. Given a 
point po £ li w e call a section of 7 at po a smooth codimension one surface S C Tie transverse 



Figure 14: Poincare map: global perturbations of a periodic orbit 7 can be studied locally 
in terms of a map ip : S — > S 1 defined by the flow of the Hamiltonian vector field Xff. 

to 7 and intersecting it at p . We then define the local map ip : S — > S* near p by letting 
p' = ^(p) be the unique point obtained by following p e S around the Hamiltonian flow 
Xh for a time close to the period T 7 of 7 (see Figured!]). Note that fixed points p = ^(p) 
(respectively periodic points p = ip k (p),k > 2) of -0 correspond to periodic orbits of the 
Hamiltonian flow Xjj of period close to T 7 (respectively close to fcT 7 ). In particular, since 
Po — ^(Po) we define the Poincare map as the differential of ip at po [2H] 



We say that the periodic orbit 7 is non-degenerate if and only if 1 is not an eigenvalue of 
the Poincare map. This is a way of saying that 7 is isolated on in the sense that there 
are no periodic orbits on Eg arbitrarily close to it. However, although 7 is isolated on Eg, 




s 
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it belongs to a continuous 1-parameter family ^ E of periodic orbits intersecting T> E at 7(see 
Figure [5]) • This is the content of the "cylinder theorem" (see for instance [2E1I2H])- 

If 7 is stable, then the eigenvalues of the Poincare map defined at a point p e 7 come 
in complex conjugate pairs of the form (e™ a , e~ tVa \ v a e M. and hence the Poincare map 
is merely a product of rotations by angles v a in n — 1 disjoint planes C T Po S. In 



Figure 15: The infinitesimal torus around a stable isolated periodic orbit 7 (p = 1) illustrated 
in the case n = 2 where there is only one stability angle v a and T Po S = R^. 

other words, every point £ 7 of the stable isolated periodic orbit 7 is surrounded by an 
infinitesimal torus Sp x . . . x Sp , where Sp = {x a e 1^ | ||x Q || 2 = F a } C M 2 , which is 
preserved by the Poincare map to first approximation in F a <C 1. By the cylinder theorem 
the periodic orbit 7 belongs to a continuous family 7^ parametrised by the energy E, and 
so one could now apply the Bohr-Sommerfeld-Maslov quantisation conditions to the family 
of torii A e 7 £ x Sp 2 x . . . x S Fn just constructed (see Figure fT5l) 



where 7 is the closed path on A consisting of a classical path going from T Po S once around 
A back to T Po S and the set of arcs of angles —v a on T Po S to close off this classical path (see 
red curve in Figure IT5|) . 

Consider the 2-dimensional surface T bounded by the periodic orbit 7 and the closed 
curve 7, constructed in the obvious way: at any point t 7^ along the curve j(t), T looks 
locally like {"f(t) + ry(t)\0 < t < T, 0<t<1} where y(t) is the transversal vector to 7 
joining the points 7(t) and j(t). At t = we complete the surface by adding the sections of 
the disc of angle —u a on T Po S. Then by Stokes's theorem we have 





a — 2, ... ,n 
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On the part of T corresponding to t 7^ we have u\r = since the tangent space to T is 
spanned by Xh and the transversal vector y {i y ix H u = i y dH = y(H) = since y lies in 
the energy surface And since r t=0 looks like sections of angle — u a of the disc of radius 
s/F^ it follows that 



a 



a=2 



On the other hand we have that 



a 



Pa 



uj = 2irF ty 



Dp 

Pat 



where D l Fa is the disc in bounded by the circle Sp a . The last equality follows by a direct 
computation, in analogy with the harmonic oscillator. Finally, by combining all the above 
we obtain [25l [26] 



a 



a=2 V 7 



h + o(h 2 ), 



(B.l) 



which is just the Bohr- Sommerf eld condition for an isolated periodic orbit. 



C Dirac brackets 

Just as in [T3], in this paper we work in conformal static gauge in order to isolate the 
physical degrees of freedom of the string. This is done by imposing the Virasoro constraints 
and static gauge fixing condition. However, these constraints together form a set of second 
class constraints and so to consistently impose these constraints from the outset one must 
work with Dirac brackets instead of Poisson brackets. In this section we show that the for 
the type of brackets (trfi(x), •} considered in section [2] this distinction does not matter since 

{trfi(x),/} D . B . = {trfi(x),/}p. B . 

for an arbitrary function / of the principal chiral model fields j = —g~ l dg and so by abuse 
of notation we drop the suffices on both brackets and write {•, •} throughout section [21 

We start with the Poisson bracket (12. 7p . To compute Poisson brackets on the circle we 
shall work on the universal cover R. So let o~\ = a + 2ir, a 2 = cr and 03 = a' in (12.71) to 
obtain the Poisson bracket {fi(cr, x)®Ji(a', x')}. This easily leads to the Poisson brackets 
{Q(a,x)fj±(a')} after noting from the definition of J±(x) that Ji(0) = — j-) and 
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lim x ^ 00 (— x) J\{x) = k(J+ m particular 



{n(a, x)®j±{a')} P . B . = (T{a + 2ir, a', x) ® 1) x - a - 2tt) - - a))^=^-^ 



47T 1 ± X 



0- + 2vr, cr) 



2tt 2x 1 

TxT^^' (x ± 1)Ji(a '' x) ® 1 ± 1 ® 2 (j+(a/) " 

x (T(a', <x, x) 



where we have used the definitions of the r, s-matrices [H] which involve the tensor product 
V — \ a a ® &a- Using the identity tr2(r/l ® A) = A for any matrix A G su(2) one can show 
that after multiplying the above equation by 1 ® j±(cr') and taking the trace tr2 over the 
second tensor factor the commutator disappears and we are left with 

1 1 47T 

fl(a, x), —trj± (cr') > = —={5{a'— a— 2ix)— 5{a' — a))T{a+2ix, a', x) J±(cr', x)T(a', a, x), 

2 J p.b. V A 

where J±(a', x) = j±(cr')/ (1 =F x). Next we multiply both sides by e ±m<J ' and integrate over 
a' from to 2ir. However, since we are on the universal cover R of S l we get two non-zero 
contributions, namely from the integrations over the two lifts [0, 2tt] and [2tt, 4tt] (assuming 
a G (0,27r)). Definition the Virasoro generators 

f 27T , line' 1 -2 / A f ^ /"^ A ' „-ina' 1 -2 /_ n 



we can write the result as follows 



5tt ./o 2^ y ' 8tt i 2~ 



{Q(a,x),L n } P . B . = -e iruT [J + (a,x),n(a,x)}, {tt(a,x), L n } P . B . = - e - m<J [J_(a, x), fi(<r, x)}. 

Note that in the above calculation it is because of the presence of the s-matrix, which arises 
from non-ultralocality of the Poisson brackets of the model, that we end up with the correct 
transformation property for Q(x) under conformal transformations. Finally, since the right 
hand sides are commutators, taking the trace shows that trQ(x) is invariant under conformal 
transformations generated by L n ,L n , namely 

{trf](x),L n }p. B . = {trfi(x),L n } RB . = 0. 

The assertion that the Dirac and Poisson brackets involving the quantity trQ(x) are equal 
now follows from the definition of the Dirac bracket which in the present case reads, 

{trfi(x), /} d .b. = {trfi(x), /}p. B . - {tr Cl(x), L n } P . B .{L n , £ m }p. B {L m , f} P . B . 

— {tr Q(x), L n } PB {L n , L m } PB {L m , /}p.b., 

for any function / of the principal chiral model fields j = —g~ 1 dg. 
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D Pinching an a-period 



In this appendix we determine the behaviour of the Riemann ^-function when the underlying 
algebraic curve E becomes singular [10, SI] • To determine the effect of degenerating an a- 
cycle on the algebraic curve E, let us consider a family E e of Riemann surfaces (e > 0) of 
genus g + 1 with homology basis {al,bl} 9 i=Q of i^E 6 ,! 
holomorphic 1-forms canonically normalised as 



Let {^iY i=0 be a dual basis of 



/ 



(d.i; 



We model the pinching of an a-cycle of the algebraic curve E by choosing a family {E e } e>0 
for which a particular marked cycle a on E e homotopic to a e shrinks to a point Pq in the 
singular limit e — > 0. The resulting surface E° is singular at Pq, and we denote by E' its 
desingularisation. 





E' 




In the limit e — > where the cycle a e shrinks to a single point Pq G E°, the cycles ao and a'o 
are homotopic to the punctures P + and P ~ on E' corresponding to the desingularisation of 

the 1-form cJq acquires simple poles 



P on E°. It follows from fID.ip that in the limit e 
at the pair of points Pq with residues 



resp+wn 

*0 



1 

2~Ki 



u 



1 



res 



p-UQ 



1 

2ttz 



Wo 



a'o 



2m 



Since u;o has no further poles it is a normalised ( f u;o = 0, « = 1, . . . , o) Abelian differential 
of the third kind on E'. Moreover, {wj}f =1 is a basis of holomorphic 1-forms on E' dual to the 
homology basis {aj,6j}f =1 for E'. Since the curve b starts and ends at P^, the component 
of the period matrix will blow up as e — > 0. All other components of the period 

f, e uj^ stay finite in the limit e — > 0. The behaviour of the 

J On J 



n 



00 



matrix Ilf 



lb, and n oj 



Riemann ^-function associated with E e 



6(z; IT) = exp <j i(m, z) + 7T2(ITm, m) 



(D.2) 



can now be analysed in the limit e — > 0. Using the fact that the imaginary part Imll of the 
period matrix II is positive definite we have Imlloo = Im (ile^, e^) > 0. It follows that 
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the quantity e mn °° is vanishingly small in the limit e — > and one finds that (ID . 2j) can be 
expanded as follows 

6(z; IP) = 6(z; n e ) + p(z + n^; W)e izo + 6{z - II*; W)e- iz °] e™ 11 ™ + O (e 2niU °°) (D.3) 
where 



z 



e c ff+1 , n £ 



il oo LL 



z J ' \ n e if 
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